JOURNAL OF RESEARCH of the National Bureau of Standards -A. Physics and Chemistry 
Vol. 77 A, No. 1, January-February 1973 

Theory of Disclinations: II. Continuous and Discrete 
Disclinations in Anisotropic Elasticity 

R. deWit 

Institute for Materials Research, National Bureau of Standards, Washington, D.C. 20234 

(August 18, 1972) 

A general theory of stationary disclinations for a linearly elastic, infinitely extended, homogeneous 
body is developed. Dislocation theory is extended in three different ways to include disclinations, i.e., 
from continuous distributions, discrete lines, and continuous distributions of infinitesimal loops. This 
leads to three independent definitions of the disclination, which can be uniquely related to each oilier. 
These interrelationships clarify Anthony and Mura's approaches to disclination theory, which at first 
appear to diverge from the present theory. Mura's "plastic distortion" and "plastic rotation" are 
identified as the dislocation and disclination loop densities. The elastic strain and bend-twist are 
derived as closed integrals in terms of the defect densities, and shown to be state quantities. The 
theory reduces to classical dislocation theory when the disclinations vanish. For every discrete dis- 
clination line, it is always possible to find a "dislocation model," which is a dislocation wall terminating 
on the line that gives exactly the same elastic strain and stress. 
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1 . Introduction 

This paper develops a general theory of stationary disclinations in a linearly elastic, infinitely 
extended, homogeneous body. 

Nabarro [l] 1 and Anthony [2] have already reviewed many of the geometrical aspects of dis- 
crete disclinations. The author [3] and Anthony [2] have also treated the geometry of continuous 
disclinations (i.e., continuous distributions of disclinations). We have a slight difference in inter- 
pretation with Anthony over the meaning of the dislocation density tensor, and the continuity 
equation for dislocations, in the presence of disclinations. This difference arises due to the essen- 
tially independent definitions of continuous and discrete defects, and we believe this problem has 
been resolved in the present paper. Furthermore, Anthony [2] has treated the statics of a wedge 
disclination in isotropic couple-stress theory, and also given a nonlinear formulation of the geometry 
of continuous disclinations. 

The development of the theory we present owes much to Mura [4]. He has developed an anisot- 
ropic theory of discrete disclinations and generalized this to continuous distributions. He does 
not try to distinguish between geometry and statics. Mura introduced the concepts of "plastic 
distortion," j3*, and "plastic rotation," <£>*, which we [3] claimed could not be defined when dis- 
clinations are present. However, these concepts can be incorporated into the general theory with 
a slight modification in interpretation: They can logically be regarded as the dislocation and 
disclination loop density tensors. Then these quantities, which have a clear definition for discrete 
defects, retain their physical significance when the theory is generalized to continuous distribu- 
tions. In this sense the dislocation loop theory of Kroupa [5] is then extended to disclinations. 

Mura also introduces an "elastic distortion" in the presence of disclinations. The physical 
significance of this quantity still has to be assessed, because it turns out not to be a state quantity, 
contrary to all previously defined elastic fields. By definition a state quantity is a continuously 
varying quantity which can be measured (in principle) by macroscopic experiments without any 
knowledge of the former states of the body. While for dislocations the elastic distortion is a state 
quantity, we find that with disclinations the elastic strain and bend-twist are the relevant state 
quantities. All of Mura's results [4] have been incorporated in the present paper. 

There are several ways in which we can divide the theory of disclinations for the purpose of 
presentation. To clarify the logical structure of this paper we discuss them next, followed by an 
outline of the paper. 

1.1. Dislocations Versus Disclinations 

We can juxtapose the two types of defects. We first discuss the known results for dislocations 
to emphasize the new results that follow when disclinations are introduced. In this sense what has 
been called "theory of disclinations" in the literature is really a combined theory of disclinations 
and dislocations. For this reason, we shall use the term defects in this paper to denote the combi- 



1 Figures in brackets indicate the literature references at the end of this paper. 
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nation of both disclinations and dislocations. 2 In this terminology, then, this paper deals with 
continuous and discrete "defects" in anisotropic elasticity. 

"Defect theory" is an extension of dislocation theory and reduces to it when the disclinations 
vanish. Although it is possible to have a pure "dislocation theory," it is not possible to have a 
"disclination theory" without dislocations. Nevertheless we left the word "disclination" in the title 
of this paper to emphasize that what is new is due to the introduction of disclinations. 

1.2 Continuous Defects Versus Discrete Defect Lines Versus Defect Loop Densities 

There are three essentially different ways to define defects: 

The geometrical theory of continuous (distributions of) defects can be formulated by examining 
the consequences of violating the classical compatibility conditions [3]. In view of this a body 
with continuous defects is also called an incompatible body, whereas a defect-free body is called 
compatible. 

On the other hand, Weingarten's theorem [3, 6] is the point of departure for the theory of dis- 
crete defect lines. This theorem provides the two characteristic constants associated with the 
discrete defect: The general Burgers vector, which reduces to the classical Burgers vector for a 
discrete dislocation, and the characteristic rotation vector of the discrete disclination. Like the 
Burgers vector for dislocations, the characteristic rotation vector plays just as important a role 
for disclinations. We have therefore ventured to call it the Frank vector, in honor of F. C. Frank, 
who coined the word "disclination." 3 

A third way to formulate defect theory is in terms of (a continuous distribution of infinitesimal) 
defect loop densities. This could be regarded as the simplest approach, because in general the 
loop densities can be arbitrarily prescribed. Furthermore any given defect can be built up from 
some loop distribution. 

In relating the three concepts in the above three paragraphs complications arise. For dislo- 
cations there is a straightforward correlation between the dislocation density, the Burgers vector, 
and the dislocation loop density. When disclinations are introduced there is an analogous corre- 
lation between the disclination density tensor, the Frank vector, and the disclination loop density. 
However, it is found that there is no unique correlation between the dislocations defined in the 
three formulations, but that they get mixed with the disclinations. This is basically the source of 
our differences with Anthony (continuous versus discrete) and Mura (lines versus loops). 

1.3. Geometry Versus Statics 

These are the stationary equivalents of kinematics versus dynamics, a distinction made in 
almost all other fields of science. 

The relations described under geometry headings simply result from the properties of Euclidean 
space and are independent of the properties of the body. The distinction between plastic and elastic 
under this heading is therefore arbitrary, but it acquires physical significance under statics. In 
every case, the results given under geometry in this paper are valid for a linearly elastic, homo- 
geneous body, finite or infinite. 

Under the statics headings the elastic quantities from geometry are related to the properties 



2 This usage would seem to ignore the possibility of point defects. However, they can he represented in this theory by discrete defect loop densities. This will be 
illustrated by an example elsewhere [24). 

3 A historical remark is appropriate here: After Weingarten |6| published his theorem, Volterra |7| recognized its implications for the discrete defects, which 
he called "distortions." This would then at first seem like the best word to use for disclinations and dislocations combined. Kroner |8| did refer to them as "Volterra 
distortions," but he subsequently used the term "distortion" for the gradient of the displacement, a usage which has now become widely accepted. Meanwhile Love 
[9] ventured to call them "dislocations." This name has stuck, but since the translational type defect played the more important role in plasticity, the name "dis- 
location" gradually became more and more associated with this type of defect. It then became necessary to distinguish between translational and rotational dis- 
locations, or, referring back to Volterra, between dislocations of the first, second, and third order, and dislocations of the fourth, fifth, and sixth order. Therefore, 
Frank [10] coined the term "disinclinations" for the latter. He subsequently modified it to "disclinations." Now the names "dislocations" and "disclinations" are 
becoming well-established in the literature, but so far no suitable term has been introduced to describe the combination of both of them, which we have simply called 
"defects." For additional information see Nabarro [3], pp. 17-20. We also add a remark on the nomenclature: In the present linear theory the rotation is traditionally 
represented by a vector, and successive rotations commute, in accordance with the usual vector addition rules. When the theory is extended to the nonlinear range 
the rotations can be finite and may no longer commute. In this case the rotation can still be represented by a vector (sometimes called rer.sor), which obeys a more 
general addition rule. See section 5.1. 
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of the body. We confine ourselves to linear elasticity, i.e., Hooke's law, and a homogeneous body. 
Furthermore the specific results in this paper are limited by boundary conditions. These are that 
the body is infinitely extended and that the displacement and strain fields vanish at infinity faster 
than r' 1 and r -2 , respectively. 

1.4. Outline of Paper 

In section 2 we derive a general solution of the plastic strain problem, which can be posed 
without specifying the nature of the defects involved. This problem is a generalization of Eshelby's 
"transformation problem" [11]. The result forms the basis for all subsequent applications to statics. 
In the derivation the general boundary conditions are introduced. Also the useful Green's tensor 
is defined. We show how Eshelby's method of solution can be generalized. Finally we show that 
a compatible plastic strain gives no elastic fields. 

In section 3 we review classical dislocation theory to set the stage for the following sections. 

In section 4 we derive the fields for a continuous distribution of defects. The geometry for this 
case has already been treated [3] and equations quoted from this reference are denoted by (II), 
(12), etc. In this section the useful incompatibility source tensor of Simmons and Bullough [12] is 
used to find the elastic strain and show that it is a state quantity. 

Section 5 treats the discrete defect line. Weingarten's theorem is used to motivate the defini- 
tion. Then the appropriate plastic strain and bend-twist are found, which are logically expressed 
in terms of Mura's /3* and </>*. The dislocation density tensor is found to depend on the Frank 
vector, an example of the mixing referred to above. The static results are derived. 

Section 6 shows the formulation in terms of a continuous distribution of infinitesimal defect 
loops. Here we identify /3* and </>*, introduced for the discrete defect, with the dislocation and 
disclination loop density tensors. Hence this gives meaning to them in Mura's generalization to 
continuous distributions. As another example of mixing we find that the dislocation density depends 
on the disclination loop density. 

In section 7 we derive some results for the discrete dipole line, and show how they are related 
to the dislocation dipole. 

Section 8 shows that the general results of sections 4-6 reduce to those of section 3 when no 
disclinations are present. 

Section 9 examines more closely the meaning of /3* introduced for a discrete defect loop. 
Without <£>*, it is shown to represent a terminating dislocation wall, i.e., a constant dislocation 
density on a surface terminating at the discrete dislocation line. This has been called the "dis- 
location model" of the discrete defect line. The elastic fields are derived. 

Section 10 examines the meaning of <f> * for a discrete loop. Without /8 *, it is shown to represent 
a compensated disclination loop, i.e., a surface with a constant dislocation density (opposite to that 
of sec. 9) terminating at a discrete disclination loop. The elastic fields of this defect vanish. 

In an appendix we develop a special notation to deal with delta functions on curves and 
surfaces. This notation is very convenient and simplifies the equations that occur when we treat 
discrete defects. 

Throughout the development of this paper we find that many concepts or quantities from dis- 
location theory generalize into pairs in defect theory. For example, dislocations generalize to defects, 
consisting of disclinations and dislocations, or the Burgers vector generalizes to the characteristic 
vectors, consisting of the Frank vector and the total Burgers vector. We have found it useful to 
introduce the new concept of "basic fields," consisting of the strain and bend-twist. Then the dis- 
tortion in dislocation theory generalizes to the basic fields in defect theory. The nomenclature that 
has developed is summarized in tables 1 and 2. 

This paper basically addresses itself to solving boundary value problems. The important sub- 
ject of the forces on and the energy of the defects introduces additional complications. It is there- 
fore omitted and will be treated elsewhere. 

Neither shall we treat applications to special problems in the present paper. These will also be 
left for future publications [25, 26]. 
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2. General Solution of the Plastic Strain Problem 
2.1. Statement of the Problem 

In this section we give a formal solution of the following problem, which can be posed without 
specifying the nature of the defects involved: given an infinitely extended homogeneous anisotropic 
body with the plastic strain ej! t given as a prescribed function of space. To find the resulting total 
displacement u' m throughout the body. 

This problem is a generalization of Eshelby's [11] "transformation problem" to an anisotropic 
medium and an inhomogeneous stress free strain. We remark here that the concept of "stress free 
strain" is identical to that of plastic strain. 

The statement of this problem can be rephrased in a manner that is almost identical to that of 
the classical problem of elasticity, i.e., in terms of the equation of equilibrium, Hooke's law, 
and the definition of strain. Using this approach we shall derive the solution to both problems 
simultaneously. 

The equation of equilibrium for the stress aij is: 4 

<rtj,i + fj = Q, (i,;=l,2,3), (2.1) 

where f\ is the body force per unit volume. Here we have used the Einstein summation convention, 
and the subscripted comma followed by the index i indicates partial differentiation: dcr^/d*;. 
The stress is related to the elastic strain e^i by Hooke's law: 

<Tij = Cijkieki, (2.2) 

where the C ijk i are the anisotropic elastic constants. Since cr,j and e k i are symmetric it follows 
that 5 

Cijki = Cjiki = Ctjue. (2.3) 



The total strain e' kl is defined by 6 



(2.4) 



In general the total deformation is not completely elastic, but part of it is stress free or plastic, so 
that 

e *i = «*'+««• (2.5) 

The above relations are conveniently combined into the expression 

CijkiuJ ki + fj = C ijk ie^ lti . 

This is the set of partial differential equations we wish to solve for uj when fj and e£, are given. 

2.2. Definition and Application of Green's Tensor 

To integrate the eq (2.6) it is useful to introduce Green's tensor function G jn (r), which repre- 
sents the displacement in the Xj direction at the field point r arising from a point force in the x n 
direction at the origin. Thus Gj n for an infinitely extended body is defined by 



4 Cartesian coordinates are used for simplicity. 

5 If an elastic strain energy function exists, we also nave the symmetry condition Ci ik i = C k uj, but we do not need this relation in the present paper. 

8 We define the symmetric part of a tensor T\ t by 7*<y)"" 1/2 (Ttj + Tjj). When T t j is a complicated expression involving many other subscripts we shall also write 
equivalently (Tij)(ij). 
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C m G jn) ik(r)+8in8(r)=0 (2.7) 

together with boundary condition that Gjn vanish at infinity. 7 Here 8i n is the Kronecker delta and 
8(r) is the three-dimensional Dirac delta function. The latter is defined in appendix B, where we 
also show that it is homogeneous of degree (—3) in r. Therefore it follows from (2.7) and the 
boundary conditions that Gj n (r) is homogeneous of degree (— 1), i.e., Gj n (r) varies as r _1 . 

We can now derive the solution of (2.6) in terms of Green's tensor. Writing for the relative 
radius vector 

R = r-r' (2.8) 

we have by (B3) and (2.7) 

= - j C ijkl G jn , ik (R)uf(r')dV, (2.9) 

where the integrations are taken over all space. Now, for any tensor T'(R), which is a function of 
R only, it follows from (2.8) that 

r,(R)=-7>(R) --ar/a*;- (2.10) 

Therefore we can also write (2.9) as follows H 

ul(r)=- j C ri ,,G jt ,j> k ,(R)u](r')dV'. (2.11) 

By the divergence theorem, (Al) in appendix A, this relation can be transformed into 

Bj(r) =-i C iJkl G jn ,r(R)^(r')dS' k + i C ijkl Gj„(R)uJ ,,(r') A' ~ f C mi G Jtt (R)uJ kli ,(r')dV , 

(2.12) 

where the surface integrals are taken at infinity. This step is usually called partial integration. 
The above relations hold only if the integrals converge. 

We now assume the following boundary condition: The total displacement itf(r) — > as r^> 

°°. Then uf k (r) will approach zero faster than r~ l as r^oo. Hence in view of the behavior of G jn , 
the integrands of the surface integrals in (2.12) will approach zero faster than (r')~ 2 as r —> o°, 
and so these integrals will vanish. With this condition we also see that the volume integrals in 
(2.9), (2.11), and (2.12) converge. Thus we have by (2.6) 

i£(r) = j G jn (R)[fj(r') -C y ,< Mf (r')]dF- (2.13) 

2.3. Solution of Classical Elasticity Problem 

In classical elasticity, which we also call compatible theory, there is no plasticity, e^ = 0, 
and therefore (2.13) reduces to 

ul(r) = JG jn (R)fj(r')dV'. (2.14) 



7 If Cijki = Ci<iij, we also have the symmetry relation Gj„ = G„j for an infinitely extended body, but we do not need this relation in the present paper. 

8 In (2,11) the Einstein summation convention also applies between the primed and unprimed indices. 
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This is the well-known classical solution, which is almost obvious if we remember the meaning of 
Green's tensor. In view of (2.6) and the boundary condition on uj, it is necessary to assume the 
following condition: The prescribed body force fj(r) must approach zero faster than r~ 2 as r— > oo. 
This requirement also insures that the integrand of (2.14) approaches zero faster than (r')~ 3 as 
r* — * °°, so that the integral is finite. 

2.4. Solution of the Plastic Strain Problem 

In the present problem, which we also call incompatible theory, there is plasticity, but no 
body force, fj = 0, and therefore (2.13) reduces to 

ul(r)=- j C m G jn (R)er fi ,(r')dV' 

= - j C Ukl Gj nJ (R)e^(r f )dV f (2.15) 

by another partial integration and (2.10). Again we have taken the surface integral to vanish in 
the partial integration, and for this it is sufficient to assume the following condition: The pre- 
scribed plastic strain ^(r) must approach zero faster than r~ l as r— > oo. This requirement is 
consistent with the boundary condition on uj in view of (2.6). Furthermore, the integrand of (2.15) 
will then approach zero faster than (r')~ 3 as r — > oo, so that the integral is finite. 

Equation (2.15) applies to any defect which can be described by the given plastic strain. It 
forms the basis for all subsequent applications to statics, where the same boundary condition on 
uZF and the condition on gf must ^^ caticfip»rl A similar rfi **olt was derived hv Mura IT31. for a time- 
dependent plastic deformation. The above derivation emphasizes that it is only necessary to know 
the plastic strain (and not the plastic distortion), and it gives the condition it has to satisfy, as well 
as the boundary condition on the displacement, to find the solution for an infinitely extended body. 

In the subsequent developments the results under geometry will hold regardless of the behavior 
at infinity, but those under statics are subject to the boundary conditions stated here because of the 
partial integrations involved. 

2.5. Eshelby's Method 

Physical science abounds with so-called "tricks" used to overcome mathematical difficulties. 
During the development of the mathematical theory of dislocations, Eshelby's bag of tricks has been 
remarkable. We wish to show how one of his recipes can be generalized to obtain (2.15). 

He calculates the displacement for his transformation problem with the help of a sequence of 
imaginary cutting, straining and welding operations [11]. He does this so that he can introduce a 
fictitious body force simultaneously with the transformed inclusion such that there is no displace- 
ment. Then, he removes the body force and finds the resulting displacement from the classical 
expression. 

To generalize this approach to our case it is presumably necessary to conceive of a continuous 
distribution of cutting, straining and welding operations. It is not important whether or not this can 
be imagined, because the formal steps are the same as in Eshelby's recipe. With the plastic strain 
e l k \ we introduce a fictitious body force given by 

fj = C ijkl e p kl .. (2.16) 

Then (2.13) shows that there is no displacement. Next, we remove this body force by applying an 
equal but opposite body force —fj= —C ijk i e p M . to the body. From (2.14) this leads to the displace- 
ment 



ul (r) = - j C m G jn ( R ) e£, ., ( r' ) dV 
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in agreement with (2.15). So the moral of Eshelby's method is that a defect described by the plastic 
strain e p k{ can be simulated by the fictitious body force given by — Cijki e P { .. 

2.6 Compatible Plastic Strain 

There is an interesting consequence of (2.15) if the plastic strain is derivable from a plastic 
displacement 

<,= <*>• (2-17) 

This will be called a compatible plastic strain. Then we have 

ul(r) = -jc ijkl G jn ,i(R) u( §k ,(r') dV 
= -[ Cij kl G jn ,i k (R) uf(r')dr 
= ]6in8(R) ufir^dF 

<(r). (2.18) 

Here the second equality follows by a partial integration, and the third by (2.7). It follows therefore 
that when the plastic strain is compatible, the elastic displacement vanishes: 

u n = u T n -u p n = 0. (2.19) 

Hence, in this case all elastic fields vanish. 

3. Review of Dislocation Theory 

3.1. Continuous Distribution of Dislocations 

3.1.1. Geometry 

If in addition to the plastic strain e£ l9 the plastic rotation (o p is also prescribed, then we can 
identify the defects as dislocations. In this case we can say that the plastic distortion (15.1) 9 

jBJ&=* + «mK (3.1) 

is prescribed arbitrarily as a function of space, where €kiq is the permutation symbol. Here e P f is 
the symmetric part of fi^ 

4i = PU (3.2) 

and the second term in (3.1) is the antisymmetric part of /3£,, or 

®J=l/2e«d8£ . (3.3) 

The dislocation density is defined by (15.2) 

a p i= — e P mkP P i iTn , (3.4) 



or, equivalently, 



9 The symbol (15.1) refers to eq (5.1) in reference [3]. 



€ P mkOCpl — p m/> k p kl m . (3.5) 
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This definition implies we are using sign convention A, or FS/RH [14|. The continuity condition for 
dislocations (15.4) 



dpi, P = 



(3.6) 



follows directly from (3.4). It implies that dislocations can not end inside the body. 

We define a Burgers circuit as any closed curve \ inside the body (fig. 1). The Burgers vector 




Defect 
Line L 



Burgers 

Circuit 
X 



Figure 1. Geometry of the discrete defect line L, its defect surface S, a Burgers circuit A., and its 

Burgers surface cr. 

The surfaces S and cr are arbitrary and curved, subject only to the condition that they terminate on L and X. The intersection of L with cr 
is labeled r /,r and X crosses 5 at r s '\ These are positive crossings by the right-hand rule. 



associated with \ is defined by Jie closed line integral 



bi=-j> ffi t dL k . 



(3.7 



This relation can be interpreted as follows: Starting with a perfect crystal we can imagine that the 
plastic deformation is produced by letting dislocations migrate into the crystal. A number of 
dislocations cut through \, producing relative displacements — /3^dL k in the lattice at the curve. 
By adding these contributions we measure the resultant Burgers vector of all the dislocations that 
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remain stuck through the surface cr bounded by A.. From (3.4) and Stokes' theorem (A2) we can also 
write 



bi=j a p ,dS p . (3.8) 



This relation shows that the dislocation density a p i represents the flux of dislocation (or Burgers 
vector) in the %\ direction that crosses unit area of a plane normal to the x p direction. 

The existence of the plastic distortion implies the existence of the elastic distortion (12.1 and 
15.5) 

Since the existence of plastic distortion implies there are no disclinations, this relation holds only if 
there are no disclinations present. This relation also allows us to draw the following conclusions 
from the definitions (3.4) and (3.7) 

£ P mk f3kl,m=a P h (3.10) 

j> f3 k idL k = b,. (3.11) 

These relations are called the basic geometric laws or field equations for a p i and bi. 
Other quantities, which we shall find useful later, are the plastic bend-twist (15.16) 

< q = < k =V 2 ««rP P r,,V (3-12) 

where the last equality follows from (3.3), the elastic strain 

€mn = P(mn), (3.13) 

and the elastic bend-twist (15.18) 

Kst = (Dt, s= 1/2 €tmn Pmn, *. (3.14) 

3.1.2. Statics 
From (3.2) and (2.3) we find that the displacement (2.15) becomes 

aS(r) =- j C ijkl G jn ,> i(R) j8£,(r') dV\ (3.15) 

where /%(r) must satisfy the condition that it approaches zero faster than r _1 as r^> oo. From this 
relation we find the total distortion 

ui, m (r)=- j C ijkl G Jn , im (R) ffi,(r') dV' 
= -jc i j kl Gj n , i (R)/3 p kl , m ,(r')dV 
= |c«wC^,,(R) U pmk a pl (r')-p p mh ,,(r')]dV 



/ 



epmkCijktGjn, i(R) ap i(r') dV'+/3^ n (r). (3.16) 
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Here the first equality follows simply by differentiating under the integral sign, where Green's 
tensor is the only function depending on r, the second equality by a partial integration, the third 
from (3.5), and the fourth by a partial integration, (2.7) and (B3). From (3.9) we obtain the elastic 
distortion 



0mn(r)=J €pmkCw£jn, i(R) Mr') dV . (3.17) 



This is Mura's half of the Mura- Willis formula f 12, 13, 15]. Willis' half would be included if (2.14) 
had also been included in the derivation. We note that from the above condition on fi^ and (3.4) 
that ctpiir) must vanish faster than r~ 2 as r— > °°. However, for the integral (3.17) to exist it is only 
necessary that ot p i(r) vanish faster than r~ l as r^> oo. 

Here we note that a p i is a state quantity because it can be measured in the present state of the 
body, (e.g., by electron microscopy or x rays). Therefore (3.17) shows that f} mn is also a state quantity, 
because it is expressed entirely in terms of a p i. On the other hand, for example, /3£ n may not be a 
state quantity, because one may not be able to measure it without knowledge of the former states 
of the body. 

3.2. The Discrete Dislocation Line 

3.2.1. Geometry 

The discrete dislocation line L is defined as the boundary of a surface S, where the material 
below S has been plastically displaced with respect to the material above S by an amount given by 
the constant Burgers vector bi (fig. 1). 

Hence, the difference between the displacement just below and above S is given by 

[u t (r)] = bi. (3.18) 

Our problem now is to find the corresponding plastic distortion. The following is a straightforward 
procedure we have developed to find it. Assume first that S is closed, enclosing the volume V . Then 

by (B7) 

u T l (r)=8(V)b l (3.19) 

represents a displacement that is equal to bi inside V and has the jump (3.18) at the surface S. The 
corresponding distortion is by (3.9) 

/B£(f)-h%(p) 

= 8, k (V)bi 

= -8 k (S)bi, (3.20) 

where we have used the divergence theorem (B24). Here S is the closed boundary of V. We see that 
the distortion is concentrated at the surface S. Since this deformation is just a rigid translation of 
part of the body, there is no elastic distortion, and therefore the above distortion is all plastic. We 
now simply generalize this expression to the open surface S of the dislocation loop: 

/8£,(r) = -MS)6/. (3.21) 

The Burgers vector for continuous distributions of dislocations was defined in (3.7). We want 
to show that this definition agrees with the constant hi introduced above. If X is any closed curve 
that encircles L in the positive sense (fig. 1), it will cross S positively at some point r vx . Therefore 
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j8£dL fc =<b8 fc (S)ML*=fc 



(3.22) 



by (3.21) and (B15). 

Now we find the dislocation density from (3.4) and Stokes' theorem (B26) 

OL p i(r) =€ p ,nk Sfr, m (S) 6/ 

= 8 p (L)bi, 



(3.23) 



where L is the closed boundary of S, i.e., the dislocation line. The vector $p(L) is the Dirac delta 
function on the curve L and it is always parallel to L. A discrete dislocation line is called screw or 
edge when the Burgers vector is parallel or normal to the line, respectively. Therefore (3.23) shows 
that the diagonal and off-diagonal components of a p i represent the screw and edge component of the 
dislocation density, respectively (see table 1). Equation (3.23) shows how to make the transition 

Table 1. Definition of various defect quantities 





Dislocation 


Name of 


component 


Quantity 
















Disclination 


Diagonal 


Off-diagonal 


Density tensor for continuous dis- 


a 


screw 


edge 


tribution of defects. 


e 


wedge 


twist 


Density tensor for continuous dis- 


y or/3* 


prismatic 


slip 


tribution of infinitesimal defect 


J or** 


twist 


wedge 


loops. 









from a continuous distribution of dislocations to a discrete dislocation line. We note that it satisfies 
the continuity condition (3.6) 



a p i jP (r) = 3, ) , p (L) 6/=0, 
by (B28). As a cross-check we also show that (3.8) is satisfied by (3.23) (fig. 1) 

I a p idSfj= I d } j(L) bidS p = bi, 



(3.24) 



(3.25) 



where we have used (B15) again. This last relation remains valid for many dislocation lines, and 
therefore can be used to show how to make the transition from discrete dislocation lines to a con- 
tinuous distribution of dislocations: For many dislocation lines the average dislocation density 
apt represents the xi component of the sum of the Burgers vectors of all the dislocation lines that 
intersect unit area of a plane normal to the x p direction. For another interpretation of a p i consider 
the result 



f a pl (r)dV=j S p (L)bidV 

= bA 3 p (L)8(V)dV 
= bij> L 8(V)dL p 



bi 



f dL p . 

JL(V) 



(3.26) 
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Here the first equality follows from (3.23), the second from (B7), the third from (Bll), and the fourth 
from (B7), where L(V) is the part of the curve L inside V only. From this expression we see that 
the average dislocation density a p \ also represents the sum of the x\ component of the Burgers 
vectors times the projected line length in the x p direction of all the dislocation lines per unit volume. 
The equivalence of the above two interpretations of the average dislocation density can also be 
shown by the methods of quantitative stereology [16]. 

We wish to point out here that there is also another type of dislocation density in wide use, 
primarily by experimentalists. This is the total dislocation line length per unit volume, usually 
designated by p: 



It is easy to show that 



\p{v)dV=\ dL. (3.27) 

Jv Jl(v) 



p(r) = t p 8 p (L) 9 (3.28) 



where t p is the unit tangent to the dislocation line. Hence the relation between the two different 
dislocation densities is 

pb 2 = t r> a I> ib, (3.29) 

from (3.23). 

3.2.2. Statics 

Now we substitute (3.21) into (3.15) to find the displacement of a discrete dislocation line. 

uT(r)=jc ijkl G jn , ,(R) 8*(S') btdV 

=J C«wCji, f i(R)6rfS; f (3.30) 

using (B12). This equation allows us to estimate the asymptotic behavior of the displacement at 
large distances from a small dislocation loop. It is simply the same as the asymptotic behavior 
of the integrand. Since Green's tensor Gjn(r) varies as r" 1 , we see that ttj(r) will vary as r -2 as 
r— >o°. Since the strain emn(r) varies as the derivative of the displacement, it will go as r~ 3 as 
r — > oo. These results are listed in table 3. A more accurate calculation of the asymptotic displace- 
ment from a small loop can also be made from (3.30) by expanding Green's tensor as a Taylor 
series in r' for a few terms and integrating over S. If this result is specialized to isotropy, we find 
the same relations as were given by Kroupa [17]. The details will be shown in a future publication 
[26]. 

We now find the total distortion from (3.30) 

<m(r) = J C ijk £ jn ,im(R) bidS f k (3.31) 

=<p e P mhCijklGjn,i(l{) bidL' p -\- I CyklGjn, i*(R) bidSm 

= j> e P mkCij k iG jn ,i(R) bidL'p+p^ir). 

Here the first equality follows by simply differentiating under the integral sign, where Green's 
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tensor is the only function depending on r, the second equality follows from Stokes' theorem (A4), 
and the third from (2.7), (B5), and (3.21). From (3.9) we then obtain the elastic distortion 

Pmn (r) =<f> epmkCijkiQjn, i (R) bidL' p (3.32) 

This relation can of course also be obtained directly from (3.17), (3.23), and (Bll). A similar result 
was derived by Mura [13] for a moving dislocation line. 

We see that the state quantity fimn can be written as a line integral along the discrete disloca- 
tion, i.e., it is expressed entirely in terms of an integral over the only regions of the body where the 
defect is localized. In general, we suggest that for a defect which is a state quantity, the necessary 
and sufficient condition for an associated field quantity to be a state quantity is that it can be written 
as an integral over the defect. For a line defect this means that the state quantity must be a line 
integral along the defect. 

Equation (3.32) also allows us to estimate the asymptotic behavior of the distortion, and hence 
the strain, at large distances from a straight dislocation line. Due to the integration it is simply the 
same as that of Green's tensor, i.e., /3mn(r) and e mn (r) vary as r _1 as r— * oo. Since the displacement 
u T n is an integral of the distortion, it will vary as In r as r^ oo. These well-known results are also 
listed in table 3. 

3.3. Continuous Distribution of Infinitesimal Dislocation Loops 

In section 3.1 the dislocation density and the Burgers vector were defined by (3.4) and (3.7) in 
terms of the given plastic distortion. A consequence of these definitions is (3.8), relating the Burgers 
vector to the dislocation density. This equation could alternatively be used to define the dislocation 
density in terms of the Burgers vector, if this quantity is prescribed in a suitable manner. It is con- 
venient to put this relation into differential form. When the dislocations are continuously distributed 
the density tensor is defined locally by 

Ab, 
apl = AS~' (S.S3) 

As we noted above, for a distribution of discrete dislocations this represents the average dislocation 
density, where Abi is the /th component of the resulting Burgers vectors of all the dislocations 
which pierce through a surface element AS P oriented normal to the x p direction at the given point. 
In a similar way Kroupa [5] introduced the dislocation loop density tensor. He defined it as 
follows: jki represents the flux of dislocation (or Burgers vector) in the xi direction that encloses a 
unit vector in the %k direction. When the loops are continuously distributed the density tensor y^i 
is a function of the position and is defined locally by 

»t-|£. (3.34) 

Again, for a distribution of discrete loops this represents the average dislocation loop density, 
where now Abi is the /th component of the resulting Burgers vectors of all the loops which are 
pierced by the line element AL k oriented in the Xk direction at the given point. 

To find the relation between the dislocation loop density and the plastic distortion, we first 
convert (3.34) to integral form 

bi = <by k idL k . (3.35) 

Now we compare this relation to the definition (3.7). The integrands can only differ by a gradient 
with respect to Xk. Therefore we can set 
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0£ =-?**+<*, (3.36) 

where uf is an arbitrary vector field. This relation identifies the plastic distortion for a continuous 
distribution of dislocation loops. With it we can find all the relations derived in section 3.1 in terms 
of loops. For example, (3.4) leads to 

OLpi = € })t nk Jkl,m- (3.37) 

This is the fundamental relation between the dislocation loop density and the corresponding dis- 
location density. Since a plastic displacement does not contribute to the elastic fields, (c.f. sec. 2.6), 
we can set u l, { — 0, without loss of generality and so we can use 

|3£=-7j« (3.38) 

for the purpose of calculating the fields of a continuous distribution of dislocation loops. This agrees 
with Kroupa's identification. Hence we can identity the plastic distortion with the dislocation loop 
density, except for a minus sign. 

We now also can give an interpretation of (3.21) in terms of infinitesimal loops. To construct 
a discrete dislocation line L, we distribute a constant density of infinitesimal dislocation loops of 
strength b\ over any surface S, whose boundary is L. This method could be taken as an alternative 
to the definition of the discrete dislocation line given at the beginning of section 3.2. So we conclude 
that for a finite dislocation loop the loop density is given by 

y k i(r)=S k (S)b u (3.39) 

where S is a surface that spans the dislocation line. The vector 8a- (S) is the Dirac delta function on 
the surface S and it is always normal to S. In a plane we have a prismatic or a slip loop according 
to whether the Burgers vector is normal or parallel to 5, respectively. Therefore (3.39) shows that 
the diagonal and off-diagonal components of y k i represent the prismatic and slip components of 
the dislocation loop density, respectively (see table 1). Equation (3.39) shows how to make the transi- 
tion from a continuous distribution of infinitesimal loops to a finite loop. By (B15) we see that it 
satisfies (3.35), which also remains valid for many finite loops, and therefore can be used for the 
transition from finite loops to a continuous distribution: For many finite dislocation loops the average 
dislocation loop density y k i represents the xi component of the sum of the Burgers vectors of all the 
loops whose surfaces are intersected by a unit vector in the x k direction. For another interpretation 
of yici consider the result 



f y k i{v)dV=\ S k (S)bidV 

= biJ8 k (S)8(V)dV 



- h k 



8(V)dS 



~ bl kn (Ek " (3.40) 

Here the first equality follows from (3.39), the second from (B7), the third from (B12), and the fourth 
from (B7), where S(V) is the part of the surface S inside V only. From this expression we see that 
the average dislocation loop density y k i also represents the sum of the xi component of the Burgers 
vectors times the projected area on a plane normal to the jca direction of all the dislocation loops 
per unit volume. The equivalence of the above two interpretations of the average dislocation loop 
density can also be shown by the methods of quantitative stereology [16]. 
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3.4. The Dislocation Dipole 

Kroupa [18] has also treated the fields of a dislocation dipole, i.e., a close pair of dislocations 
with opposite Burgers vectors. In this section we wish to present some general formulas for such a 
defect. 

3.4.1. General Definition of a Dipole 

We first give a very general definition of the dipole conjugate to any defect, following a similar 
line of reasoning as Kroupa used. The dipole is composed of two parts: the first is obtained from the 
basic defect by translating it through a small distance, and the second is the negative of the basic 
defect at its original position. To give this concept a mathematical formulation, let the basic defect 
be given as a function of position r by the source function S(r). For example, this could be e 1 ^ 
of section 2, or a p i of section 3.1. If this defect is rigidly translated through a distance £, the source 
function of the new defect configuration becomes S(r — £). Hence the source function of the 
conjugate dipole is given by 

S"(r)=S(r- £)-S(r). (3.41) 

This result applies to a finite dipole. 

For a discrete defect concentrated on a point, line, or surface, it is customary to deal with the 
infinitesimal dipole. It is obtained by letting £ approach zero and the strength of the basic defect 
approach infinity in such a way that the field of the dipole remains finite. If £is infinitesimal, we can 
use Taylor's expansion to write (3.41) as follows: 

S D (r)=-£iSAr). (3.42) 

We see that our definition has the opposite sign from Kroupa's, but it agrees with the convention in 
electrodynamics. We shall now show that a similar relation holds between the fields of a defect and 
its conjugate dipole. Let the field of the basic defect be given by the generic expression: 

/(r)=|c(R)S(r')rfr, (3.43) 

where G(R) is some kernel of integration. For example, (2.15) and (3.17) have this form. The corre- 
sponding field of the conjugate dipole is 



f"(r) = JG(R)S'>(r')dV' 

= -tiJG(R)S,v(r')dV 
= -{ i JG, i CR)S(r')dV 



=-&/;,(r). (3.44) 

Here the second equality follows from (3.42), the third by a partial integration, and the fourth from 
(3.43). 

Equations (3.42) and (3.44) are the fundamental relations between functions of the basic defect 
and the corresponding function of its conjugate dipole. The basic defect here is arbitrary. It could 
for example be a dipole itself; in this case we obtain the dipole of a dipole, or a quadrupole. In this 
way all higher order multipoles are defined. 
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3.4.2. Application to the Discrete Dislocation Dipole 

We now apply the above results to the case of a discrete dislocation line, figure 2. For example, 
the basic equation for the dislocation density is given by (3.23). Therefore the dislocation density 




FIGURE 2. The discrete defect dipole line conjugate to the basic defect line of figure I. 



inis figure snows me nnne case, r or tne innnuesimai case £— *u. 



of the conjugate dislocation dipole is by (3.42): 

a pi (r) = ~"£nttp/,m(r) 



(3.45) 



With this explicit expression we can clarify the meaning of the infinitesimal displacement f m : We 
leg f m —* and b\ — » °° in such a way that big „, remains constant. 

Next we find the displacement field of a dislocation dipole from (3.44) 



<(r)= -£/»<„» 



— — Q) €pmkCijklGjn,i(R) bi£ m dL p + S m (S) bn£m, 



(3.46) 



where we have used (3.31) and (3.21). It is interesting to note that the last term in (3.46) can be made 
to vanish by choosing S to lie along £, or t; m dSm = 0, i.e., by letting S be an infinitesimal strip con- 
necting the basic dislocation to its displaced partner. As a consequence we see that the displace- 
ment of a dislocation dipole is a state quantity. 

From (3.46) other static quantities can be derived, such as the distortion ffi nn . The same result 
can alternatively be obtained by substituting (3.45) into (3.17). 

Equation (3.46) also allows us to estimate the asymptotic behavior of the displacement at large 
distances from a straight dislocation dipole line. By comparison with (3.32) it is the same as that 
of the distortion of a straight dislocation line, i.e., it varies as r _1 as r— » °°. Hence the strain will 
vary as r~ 2 as r— » °°. These results are also listed in table 3. We note from the table that the disloca- 
tion dipole can be classified as a defect with properties between those of the dislocation line and the 
dislocation loop. 
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Sometimes a narrow elongated dislocation loop is regarded as a dipole. This is not possible 
with our definition, since the two components of our dipole must be disjoint. However, two parallel, 
infinite, straight dislocations with opposite Burgers vectors are included in our definition of a dipole. 
This dipole resembles an elongated loop. By contrast, it is possible to regard a dislocation dipole 
as a special type of loop. 

4. Continuous Distribution of Defects 

As discussed in section 1.1, we shall denote the combination of dislocations and disclinations 
by the word defects, see table 2. 

TABLE 2. Generalization from dislocation theory to defect theory 



Quantity 


Dislocation theory 


Defect theory 




Dislocation density a 
Plastic distortion fi 1 ' 
Burgers vector b 
Elastic distortion p 
Displacement jump [u] 
Dislocation loop density yifi 1 ') 


Dislocation density a 
Disclination density 

Plastic strain e'' 




Characteristic vectors 


Plastic bend-twist k? 
Total Burgers vector B 




Frank vector ft 
Elastic strain e 




Elastic bend-twist k 

Displacement jump [u] 
Rotation jump [o>] 

Dislocation loop density y(/3*) 
Disclination loop density £ (</>*) 


Defect loop density tensors 



4.1. Geometry 

Consider now an infinitely extended body in which the plastic strain e% t and bend-twist k 1 ' are 
given as independently prescribed functions of space. For convenience we refer to these two func- 
tions as the basic plastic fields. Then we may have disclinations as well as dislocations. The dis- 
location and disclination densities are denned by (16.3, 16.1) 



a pl — ~ Cpmk (ejfci, m + tklqK^q) , 

Opq — ~ €pmkK kq ^ w , \6qp— ~ e qnlK [pi n ) • 

The continuity equations (16.8, 16.2) 

Ctpl,p+€i PQ Opq= 0, 

Upq, p ^ 



(4.1) 
(4.2) 

(4.3) 

(4.4) 



follow immediately from (4.1) and (4.2). The former means that dislocations can only end on dis- 
clinations, and conversely, if the disclination density is asymmetric, dislocations must emerge from 
it. The latter shows that disclinations cannot end inside the body. 

Associated with the Burgers circuit X we now define the characteristic vectors, the general 
Burgers vector Bi and the characteristic rotation vector £l q , which we have agreed to call the Frank 
vector in section 1.2, by 
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B, = -j> (<•[',- e„ ; ,< ( .r r )dL k , (4.5) 

n, = -^K^L fc . (4.6) 

These relations can be interpreted as follows: Starting with a perfect crystal we can imagine that 
the plastic deformation is produced by letting defects migrate into the crystal. A number of them cut 
through X. Every disclination that cuts through X produces a relative rotation — k 1 ' dLk in the lattice 
at the curve. These contributions added around the contour X give the resultant Frank vector of 
all disclinations that remain stuck through the surface cr bounded by X. The relation (4.5) is easily 
seen to be identical to (3.7), when k 1 ' is given by (3.12), by doing a partial integration, where the 
integrated part vanishes, and using (3.1) or (3.2). Therefore, it is simply a generalization of (3.7). 
It represents the resultant Burgers vector of all defects (dislocations and disclinations) that have 
cut through X and remain stuck through cr. We wish to emphasize that Bi does not represent the 
relative displacements added around the contour X. The latter quantity is not an invariant for the 
curve X, but depends on the point where the integration is started. The definition (4.5) is essentially 
motivated by Weingarten's theorem (section 5.1). By Stokes' theorem (A2) we have from (4.5-6), 
(17.4,17.3) 

Bl = — I € pJ nk (€% lt m ~ £lqmK kQ ~ €lqrK kq ^ m X r ) ubp 

= \ (otpi — eiqrOpqXr) dSp, (4.7) 

ilq=— I €pmkK kqm dbp 

J(T 



— I dpqdSp, 



(4.8) 



using (4.1) and (4.2). Relation (4.7) shows that the disclination density 6 pq also contributes to the gen- 
eral Burgers vector, in addition to the dislocation density ol p i. Relation (4.8) shows that the disclina- 
tion density 6 pq represents the flux of disclination (or Frank vector) in the x q direction that crosses 
unit area of a plane normal to the x p direction. 

For later application we also write (4.1-2) in their equivalent forms 

tpmk&pl = ^ml, k e kl, m ' €pmkKi p €klmK •> (4.") 

€pmkOpq = K mq k ~ K kQf m , \£qnld qp ~ K np / — K\p,n). (4.10) 

The incompatibility tensor is defined by (14.1, 16.6) 

7)pq — €pmk£qnie klf mn (4. 1 1 ) 

= — (€qnlOLpl, n + 6pq)(pq), (4.12) 

where the second equality follows from (4.1~2). The continuity equation for the compatibility (14.2), 

Vpq,p=0, (4.13) 

follows immediately from (4.11). 
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With disclinations we claim that the total distortion is no longer the sum of an elastic and 
plastic part, (3.9), simply because /3^ n is not denned. For, if j8^ n existed, the plastic bend-twist 
would be the gradient of the plastic rotation, (3.12), and consequently the disclination density 
(4.1) would vanish. 10 Instead, since the plastic strain e kl and bend-twist k^ q are prescribed, we 
postulate the existence of the elastic strain e/ci and bend-twist Kkq, which are called the basic elastic 
fields. So we have (12.3, 14.3, 12.16, and 15.17): 



c T = 



0) 



U {n,m)~ e mn + e mn , (4.14) 

0j fS = 1/2 e tmn ul tms = K 8t + K P st . (4J5) 

Since it follows from (4.14-15) that 

€pmk (ell, m + *klq Kl lQ ) =€ pmk U T Lkm =Q, (4.16) 

these relations allow us to derive the following basic geometric laws or field equations for a p i 
and pq from (4.1-2), (16.11, 16.10) 

€pmk (ekl,m+ €klqKmq) —Ctph (4.17) 

€pm k Kkq,m = Opq. (4.18) 

The geometric meaning of these equations is that the defects (dislocations and disclinations) 
are the sources of the basic elastic fields (elastic strain and bend-twist). It also follows from (4.14- 
15) that 

e^ — eiqr K T kq x r = (uj—u[ lr] x r ) ,k , (4.19) 

and hence these relations also allow us to derive the basic geometric laws for B\ and il Q from 
(4.5-6), (17.2, 17.1) 

<V {ekl—eiqrKkqXr) dL k =Bi, (4.20) 

(J) KkqdL k = ftq. (4.21) 

'A 

These may be regarded as the field equations in integral form, equivalent to (4.17-18). Finally, 
the relation (4.14) allows us to derive the basic geometric law or field equation for the incompati- 
bility 7] pq from (4.11), (1 4.4) 

— €p?nk€qniekl,mn = 7]pq, (4.22) 

which shows that we can isolate the incompatibility as the source of the elastic strain. 

4.2. Statics 

This section extends the work of reference [3], which gave only the geometry of a continuous 
distribution of defects, reviewed in section 4.1. The main result we shall find is an explicit and new 
expression for the elastic strain, eq (4.29), in terms of defect densities. We shall also find the elastic 
bend-twist. 



10 This situation is analogous to the case of pure dislocations: There the total displacement is no longer the sum of an elastic and plastic part, u[=u/ + «f. For, 
if uJf existed, the plastic distortion would be a gradient, /3jp,= uf k , and the dislocation density (3.4) would vanish. 
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We find the total distortion from (2.15) 



ulJr)=- j C ijkl G jnJt AR)e^,(r')d^ 
= -jc Ukl G jt ,j(R)e>:. Un ,(r')dV' 

= I C ijk iG jn ,i(R)[€ P mkOi p i(r) — e£j ifr ,(r') -ep m kKf p (r f )]dV' 
= fepnucCmiGjnAmccpiir^-KfpirWV' + e^ir). 



(4.23) 



Here the second equality follows by a partial integration, the third from (4.9) and (2.3), and the fourth 
by a partial integration, (2.7) and (B3). We see that for the integral to converge it is only necessary 
that api(r) and /c£(r) vanish faster than r _1 as r^ oo. From this relation we now proceed to derive 
the basic elastic fields and show that they are state quantities. 
The elastic strain is found from (4.14) 11 



e m n(r)= I e P mkC ij k iGjn,i(R)[a p i(r' ) — K^ p {r')\dV\ mn ). 



(4.24) 



To show this is a state quantity we shall use the concept of an incompatibility source tensor, intro- 
duced by Simmons and Bullough [12]. In contrast to Simmons and Bullough, who derived several 
forms of it from a general definition, we define it directly as follows 



Imnp<i(r) = (477) 1 J 6pmk€</slC ijklGj n j>s>(r f )R x dV\ mn ). 



(4.25) 



We see from this relation that the homogeneity of Imnpq(r) is the same as that of Gj„ (r), i.e., of 
degree (—1), so that they both vary as r _1 . By using the identity 



€pmk€qsl : 



Opq Ops Opt 

Omq Oms Oml 

$kq $ks 8kl 



(4.26) 



and (2.7) it can also be written as 

Imnpq(r) = (47r) '^(drnnSpq — &mq$np)r~ l ~ (C ijpq — C ijkl^pq) I Gj n , i'm' {r')R~ x dV 

+ (Cijmq — Cijkk$mq) I Gj n j> p ' (r' )R ~ l dV'](mn)' (4.27) 



Except for the last line, this expression agrees with the one Simmons and Bullough have called the 
Eshelby-Eddington formula. Now the incompatibility source tensor satisfies the following relation- 
ship: 

£qsllmnpq,s(l t )= (47t) _1 I € pm k[C ijklGj,, j' 8 > (r ) — Cjjk 8 Gj„ , /'/' (r' ) ]R ~ s l dV' { mn) 

= (47T)- 1 I €pmdCij k iGjn^(r f )R^l — Cij k8 Gjn t v8'(r , )R- l 1 ]dVl mn) 

= -jep m k[C ijk iG jn ^(r f )8(R)-8kMr , )^7rR);l]dVl mn) 

= — [€pmkCijklGj u ,i(r)](mn)' (4.28) 



1 See footnote 6 on page 53. 
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Here the second equality follows by partial integrations, and the third from R~* s = — 4>7t8(R) and 
(2.7). From the above the term in (4.24) containing k{ p becomes 

I eqslImnpq^i^K^ir^dV = I €q 8 Jmnpq(ll) Kf p>s ,(r')dV' 

= -jlmnpq( 1 R)e qp (r f )dV f 

by a partial integration and (4.2). Note that in the partial integration the surface integral vanishes 
because of the asymptotic behavior o{I mnP q and Kf p , discussed above. Hence we find for (4.24) 

e mn (r)= I e pmk C ijk iG jn j(R)a p i(r')dV[ mn) - I I m npq(R)Oq P (r')dV'. ^ 2 9) 

This is the result we seek. We see it is a state quantity because it depends entirely on the defect 
quantities a p i and 6 QP . 

The incompatibility source tensor Imnpq was originally introduced by Simmons and Bullough 
to solve the so-called incompatibility problem, i.e., to find the elastic strain e mn when the incom- 
patibility 7) PQ is given as a prescribed function of space. We show here how this can be done. From 
(4.28) 

Gprk^qsl' mnpq,rs\Y) \,^ijkr^jn t im\V) ^km^ijrl^jn, ir (**) J(mn) 

= ~ [CijklOjn,im(r) +dkm&ln&(r)](mn), (4.30) 

using (2.7). Therefore we find 

emn(r) = — I C lj A/G jn ,/ w (R)e^(r')^F ( ' m „-e;; // (r) 

= I [e P rk€q S iImnpq,rs(R) + d km 8i n 8(R)] e£,(r' ) dV — e£ n (r) 
= I eprkeqsiJmn P q(R)e k , lr , s ,(r')dV' 

= \Imnpq(R)Vpq(r')dV f . (4.31) 

Here the first equality follows from (4.14) and (2.15), the second from (4.30), the third from a partial 
integration, and the fourth from (4.11). 

We next wish to derive the elastic bend-twist. First we find the derivative of the total distortion 
(4.23) 

<m>)=/w^^ (4.32) 

The Ki p term in this expression becomes by a partial integration 

- I epmkCijkiGjnA^Klp^ir^dV' 

= j eprnkCijkiGjnARneqsiOqpir^-K^ir^W' 

= I €pmkCqslCijklGj n ,i(l{)6qp(r')dV' + € P mnK^,(r). 

Here the first equality follows from (4.10), and the second by a partial integration and (2.7). Hence 

u l, ms( r ) = I ZpmkCijkiGjnAsiJVlOLpiiY^dV' + I 6pm k eq S iCij k iG jn ,i(R)0qp(r ')dV 

+ <n, s (r) + € pmw < p (r). (4.33) 
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Finally we have from (4.15) 

K s t(r) = ll2 I €tmn€pmkCijklGjn, is (R ) OL p i ( r' ) dV + 1/2 I €tmn€pmk£qslC ijklGjn,i(R)6qp{v' )dV' , 

(4.34) 

which is seen to be a state quantity too. 

This section then has extended the results for dislocation theory of section 3.1.2 to a more 
general defect theory. The central result there, eq (3.17), has been replaced by (4.29) and (4.34). 
We note that concepts, quantities, or equations from dislocation theory generalize into pairs of 
concepts, quantities, or equations in defect theory. Some of these ideas have been summarized in 
table 2. 

5. The Discrete Defect Line 

5.1. Wei ngarten's Theorem 

The point of departure for the discrete defect line is the following theorem [3]: 
WEINGARTEN's THEOREM: On following around an irreducible circuit in a multiply-connected 
body, the displacement and rotation change by an amount that represents a rigid body motion, if 
and only if the classical elastic compatibility conditions are satisfied throughout the body. 
Explicitly these changes are given by (13.4, 13.3). 

[ui] = Bi-\-€i qr rt Q X r , (5.1) 

[o>J=n 9 , (5.2) 

where the constants Bi and Cl Q are given by line integrals along the irreducible circuit A. (13.6, 
13.5): 

Bl = j> (ekl — €lqrKkqX r )dL k , (5.3) 



£la = 



<j> K kq dL k . (5.4) 



That these quantities are constant is easily shown by Stokes' theorem and the compatibility 
equations, i.e. (4.17-18) with a p i = 0, 6 pq =0. We note incidentally that the definitions (5.3-4) are 
consistent with the relations (4.20-2 1). 12 

5.2. Geometry 

The discrete defect line L is defined as the boundary of a surface S, where the material below 
S has been plastically displaced with respect to the material above S by an amount which represents 
a rigid motion (fig. 1). 

Hence, the difference between the displacement just below and above S is given by 

[ui(r)]= bi+ €t<pn q (xr-x° r ) , (5.5) 

where bi represents a rigid translation and the second term a superposed constant rotation of 
amount Q, q around an axis through x Q r The constant bi will be called the Burgers vector for the 
discrete dislocation line contained in the defect line, and is to be distinguished from the general 



12 For the nonlinear generalization of the present theory it would be necessary to determine if Weingarten's theorem still holds. Then for a finite rotation (5.1) would 
have to be modified into [1/ /] = /?/+ (*iqr£lqX r + iP-flqXq — £ft 2 xi) (1 + iH 2 ) -l . Here the Frank vector or versor has the direction of the rotation axis and the magnitude 
2tan±<J>, where <l> is the angle of rotation. The addition rule becomesn,,^ (ft (l) + n ( fl 2) — teiqrCl^W, 2 *) (1 — in (,) n ( f 2) ) -' for a rotation (1) followed by a rotation (2). 
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Burgers vector defined by (4.5). The constant Cl q will be identified with the Frank vector (4.6). The 
relation (5.5) implies there is also a jump in rotation across S given by 

[o) Q (r)]= l/2€kiq[u>i(r )],k = fl 9 . (5.6) 

If we visualize the material in a tube around L removed (fig. 1), we have a doubly connected 
body, and since the jump across S represents a rigid body motion, Weingarten's theorem applies. 
Hence, this doubly connected body is compatible, i.e., the basic elastic fields satisfy the compati- 
bility equations, even on S. Therefore, the results (5.5-6) should also follow from (5.1-4). The 
real, incompatible, simply connected body is obtained by letting the cross section of the tube 
vanish. Then (5.3~4) become identical to (4.20-21), and by the compatibility of the total deformation 
these relations are equivalent to (4.5-6). 

Our problem now is how to embody the statements (5.5-6) into definitions for the basic plastic 
fields, i.e., the plastic strain and bend-twist. We give a straightforward operational procedure to 
obtain these quantities and then verify that they are correct by (5.1—2) and (4.5-6). Assume first 
that S is closed, enclosing the volume V. Then by (B7) 

uf(r) = 8(V){b / + e, q? n q (xr-x r )} (5.7) 

represents a displacement, which is the same as (5.5) inside V and vanishes outside V. Thus, it 
has the required jump across S. Equation (5.7) could be regarded as describing a grain of volume 
V and boundary S, whose orientation with respect to the surrounding material is given by the rigid 
motion (5.5). We assume that the deformation (5.7) comes about by means of a plastic deformation 
on the surface S and an elastic translation and rotation in V. To find the basic plastic fields we just 
calculate the basic total fields, i.e., the total strain and bend-twist. In general these can be split 
into elastic and plastic parts, (4.14-15). But since the only elastic deformation is a rigid motion, 
the basic elastic fields will vanish and the plastic fields will equal the total fields. We have from 
(3.9) and (5.7) 

Pl l (r)=dAV){b, + e, qr n q (xr-x°)} + d(V)ei qk (l q 

= — 8k(S){bi-\-€l Q filq(Xr — ^)}-\-8{V)€kl(flq, (5.8) 

using the divergence theorem (B24). From this we find 

ej^(r) =- dk(S) {bi + ei qr £l q (x r - x°)}(kih (5.9) 

wJ(r)=-l/2€*/,8*(S){6i+ ctvtlq(x r -aP r )} + 8(V)a q . (5.10) 

We see that (5.10) represents a rotation, which is the same inside V as (5.6) and vanishes outside 
V. Thus, it has the required jump across S. We next find the bend-twist from (4.15) 

K'j nq (r)=— 1/2 6kl Q [Sk(S){bi+ eiqr£lq(x r — X*)}], m — 8m(S)Slq, (5J1) 

using the divergence theorem again. We see that the basic total fields are concentrated at the 
surface S. As mentioned before, since the deformation is just a rigid motion of part of the body, 
there are no basic elastic fields, and therefore the fields are all plastic. The next step is to assume 
that (5.9) and (5.11) hold in the same form for the open surface S of the defect loop. To write down 
the final results it is convenient to introduce the "plastic distortion" and "plastic rotation" defined 
by Mura [4]: 

P*,(r) =-8 k (S){h+€tvn q (Xr-X*)}, (5J2) 

***M =-8k(S)£l q . (5.13) 
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These quantities will be interpreted in section 6 as the dislocation and disclination loop densities. 
Until then they serve as convenient intermediate quantities for the purpose of calculation. Now 
we write 

ei:,= /3*, , (5.14) 

<, = l/2 € ,,^*,„ +<»*„, 

(Kf,= l/2€«^8* fi , + ^). (5.15) 

These are the results we looked for. Equations (5.14-15) together with (5.12-13) represent the basic 
plastic fields for a discrete defect line L spanned by the surface S. From (5.14-15) we also have the 
result 

<,,„, + «'■"'<„ = fill . m + W^ (5.16) 

which will be useful for later purposes. As we mentioned above, we now check the validity of 
results (5.14-15). From (4.5) 



= <p 8 k (S){bi-ei qr ft q x r }dL k 



= bi-ei qr £lqX% (5.17) 

by a partial integration, (5.12-13), and (B15). This then is the relation between the general Burgers 
vector and the dislocation Burgers vector. Next we have by (5.15), (5.13), and (B15) 



-<j) K* q dL k = -j> <t£ q dL k 

= <j>8 k (S)agdL k 



= 0>q. (5.18) 

From (4.6) this relation identifies ft 9 as the Frank vector. If we substitute these results into (5.1-2) 
we obtain (5.5 -6), as was required. 

We now find the dislocation density from (4.1) and (5.16) 

oc v i{v)=-e pmk {^ lm ^e klq ^l q ) 

= epm k [5 k ,m(S){bi+ei qr O<<i(Xr — Xr)} + d k (S)€l qr Q q 8rm+€ k lq8rn(S)rL q ] 

= 8 p (L){b,+ €i qr n q (xr-x° r )}, (5.19) 
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where the first equality follows from (5.16), the second from (5.12-13), and the third by Stokes' 
theorem (B26) and a cancellation. The disclination density is obtained from (4.2) 

upq\Y) = €pmk<p kq, m 

= €pmkuk, m\b)lLq 

= 8 p (L)a q (5.20) 

where we have used (5.15), (5.13), and (B26). Relations (5.19-20) represent the defect densities for a 
discrete defect line. We see from (5.19) that the Frank vector il Q also contributes to the dislocation 
density, in addition to the dislocation Burgers vector bi. In these relations L is the closed boundary 
of S. The vector 8 P (L) is the Dirac delta function on the curve L and it is always parallel to L. A 
discrete disclination line is called wedge or twist when the Frank vector is parallel or normal to the 
line, respectively. Therefore, (5.20) shows that the diagonal and off-diagonal components of 6 pq 
represent the wedge and twist components of the disclination density, respectively (see table 1). 
Equations (5.19-20) show how to make the transition from a continuous distribution of defects to a 
discrete defect line. 

The discrete defect line, characterized by the dislocation Burgers vector bi and the Frank 
vector £l q , was defined in this section independent of the defect density tensors a p i and d pq defined 
by (4.1-2) in section 4.1. The question arises whether there is a unique correlation between these 
definitions. For dislocations only there is a straightforward relation between the Burgers vector and 
the dislocation density of a continuous distribution, given by (3.8), or between the dislocation density 
for a discrete line and its Burgers vector, given by (3.23). When disclinations are introduced there is 
a similar straightforward relation between the Frank vector and the disclination density of a con- 
tinuous distribution, given by (4.8), or between the disclination density for a discrete line and its 
Frank vector, given by (5.20). However, as we noted, the general Burgers vector now contains a 
contribution from the disclination density, eq (4.7), or the dislocation density for a discrete defect 
line contains a contribution from the Frank vector, eq (5.19). What this means is that the definitions 
of the "dislocation" are not identical in both approaches. Therefore, as we see from the relations 
quoted, in the transition from one formulation to the other a certain amount of mixing occurs. 
Anthony [2] has handled this difficulty by adopting the discrete disclination loop as the true definition 
of the disclination. Then our dislocation density tensor a p i in (5.19) partly describes the discrete 
disclination. He therefore divides a p i into two parts, a true component that corresponds to the dislo- 
cation line with Burgers vector bi, and a component that belongs to the disclination line with Frank 
vector fl Q . Hence our difference with Anthony merely reduces to a difference in point of view. 

We prefer to retain our own point of view with different definitions of the "dislocation" be- 
cause it is in fact difficult to identify the dislocation line in a discrete defect line, as will be seen 
from (5.19). For example, let us change the axis of rotation and the dislocation Burgers vector to 

< = *•+&, (5.21) 

b'^bi + eup-tlqtr. (5.22) 

Then a p i in (5.19) is unaltered. So the dislocation Burgers vector bi is not uniquely defined, but de- 
pends on the location of the axis. On the other hand, the general Burgers vector Bi in (5.17) is un- 
altered by (5.21-22), and, therefore this is the quantity which is invariant for a discrete defect 
loop. This is another motivation for introducing it. 

The significance of (5.19) can further be illustrated as follows. Consider an infinitesimal 
volume dV centered about some point on L. The jump across S is given by (5.5). The rotational 
part of (5.5) can be approximated by a constant inside dV, since it is so small, i.e., locally we can- 
not determine if the displacement jump is due to a rotation or a translation, even if we know the 
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rotation jump locally. Now the dislocation density is a local tensor field. Therefore the dislocation 
density tensor can be found by replacing hi in (3.23) by the jump (5.5), in agreement with (5.19). 
Hence, the dislocation density (5.19) at a point on the line L is exactly what would be expected on 
the basis of the local plastic displacement near the point, or a small Burgers circuit around L at 
the point. We shall illustrate this point more explicitly by examples in future publications [25, 26]. 
From (5.21-22) it would seem at first that we could eliminate the discrete dislocation from 
(5.19), by choosing f r such that bi=0, but this is only possible if b\ is normal to Cl Q . However, we 
can draw the following important conclusion from (5.19). The axis of a discrete disclination line 
can be translated from the point x° r to the point x%' by adding a discrete dislocation to the line 
with a Burgers vector given by 



bi= ei qr VLq{x r — X° r 



(5.23) 



which is normal to the Frank vector Q Q . In other words, we can move the axis by adding the dis- 
location density 

a p i(r)=8 p (L)e lqr £l Q (x o r -x r '), (5.24) 

as is evident from (5.19). 

The above development suggests the possibility that a dislocation could end on or originate 
from a disclination line. Consider the three curves L, //, and L", illustrated in figure 3, which 




Figure 3. Defect lines which join at nodes. 

meet at nodes. These curves could represent the discrete defect lines described by the following 
expressions: 

a p /(r) = 5 p (I){^ + e/ 9r n g U r -^)} + 8 p (L'){^-hez 9r a (? U r -^')} 



Mr)=S p (L + L')0,. 



d P (L")e lqr n Q (x r -x° r '), (5.25) 

(5.26) 



There are several ways in which one can view this defect. First, it consists of a defect line as 
described by (5.19-20) along the curve L and L' to which a dislocation line with Burgers vector 
(5.23) has been added along L' and L". Second, it is a dislocation line along L" , which connects 
two points of a defect line along L and L '; as a consequence the axis goes thru x® for L and thru 
x Q r ' for V. Third, it consists of two defect loops with the same Frank vector, one along L' and L" 
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with axis thru x% , and one along L and minus L" with axis thru x° r \ along the line of overlap L" the 
rotational parts cancel and only a dislocation line is left. So we see how the defect described by 
the relations (5.25-26) can be regarded as a combination of two simpler defects. We also see here 
how a discrete dislocation line can end on a discrete disclination line. A special case of this geom- 
etry has already been discussed [19] and a detailed mathematical analysis of this particular example 
will be presented in a future publication [25]. 

We note that (5.19-20) satisfy the continuity equations (4.3-4) by (B28): 

Otpl,p+ €lpqO PQ = 8p(L)€lqr£lq8rp+ €lpq$p(L)n q = 0, (5.27) 

pq ,p=8p, p (L)n q =0 9 (5.28) 

Similarly, it is easy to show that (5.25-26) also satisfy the continuity equations, confirming an 
assertion we have made [19]. 

As a cross-check we also show that (5.19-20) give consistent results for (4.7-8) 

Bi=j 8p(L){bi-€u^rf}dS p = bi--€i< p n& r , (5.29) 

I pq dS P = J dpmtlqdSp = £l Q , (5.30) 

Ja J <r 

by (B15), in agreement with (5.17-18). These relations remain valid for many defect lines as well, 
and could therefore be used to show how to make the transition from discrete lines to a continuous 
distribution of defects. However, the following relations are more convenient to make this point: 



I otpidSp=\ bp(L){bi~\- ei qr ftq(xr — x° r )}dL p 



L 



= bl — eiqrflqix 1 ^ ~ X° r ) 

= [«f] Ur , (5.31) 

0pqdSp=[o) q y (5.32) 



Here the first equality follows from (5.19), the second from (B16) where x L r a is the point of inter- 
section of the curve L with the surface cr (fig. 1), the third from (5.5), and the last from (5.6). So 
for many defect lines the average dislocation density a p i represents the xi component of the sum of 
the displacement jump vectors and the average disclination density 6 pq the x q component of the 
sum of the rotation jump (or Frank) vectors of all the defect lines that intersect unit area of a 
plane normal to the x p direction. 

5.3. Statics 

In this section we find the basic elastic fields, i.e., the elastic strain and bend-twist, for a dis- 
crete defect line. 

If we substitute (5.14) and (5.12) into (2.15) we find the displacement 

<(r) = |c^G jw , i (R)8 fe (S'){6 / +€ ig rn 9 (<-^)}dr 

= j s C m G jn AK){bi+€iMxl-4)}dS' k (5.33) 
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by (B12). This is the expression for the total displacement due to a finite discrete defect loop. It 
allows us to estimate the asymptotic behavior of the displacement at large distances from a small 
disclination loop, bi = 0. Since Green's tensor Gj n (r) varies as r _1 , we see that uj x {r) will in general 
vary as r~ 2 as r^> oo. However, for special cases such as a symmetric loop centered on its axis, the 
integral vanishes to first order by symmetry. Therefore a finite symmetric disclination loop has 
the anomalous asymptotic dependence of r~ 3 as r— ► » Since the strain e mn (r) is obtained from 
the derivative of the displacement, it will vary as r~ 4 as r— » » These results are listed in table 3. 

Table 3 

Asymptotic behavior of the displacement u T and the strain e at large distances from certain defect configurations. 
The numbers in parentheses refer to the equations from which the estimate is obtained. A symmetric disclination loop is 



one for which (x' r ~ 


-*><»;- 


0. 














Dislocation 


Disclination 




u T 


e 


u T 


e 


Line 


lnr 

r- 1 (3.46) 

r- 2 (3.30) 


r l (3.32) 

r -3 


rln r 

lnr (7.6) 

General: r -2 (5.33) 
Symmetric: r -3 (5.33) 


lnr 

r- 1 

r- 3 
r" 4 


(5.37) 


Dipole 




Loop 









For infinitesimal symmetric disclination loops both uj t and e mn vanish. A more accurate calculation 
of the asymptotic displacement from a small but finite disclination loop can also be made from 
(5.33) by expanding the Green's tensor as a Taylor series in r' for a few terms and integrating 
over S. The details will be worked out in a subsequent publication [26], where it will be shown 
that for the isotropic case the results reduce to those of Li and Gilman [20]. 
We now find the total distortion from (5.33) 

<,nM)=j^C ij uG jn j m (R){b l ^e lQr n (i (xl-x t )}dS' k 

= ie pmk Ci jk £j n j(R){bi+€i qr n q (xl.-x r )}dL^ 

+ I CijkiGjn f i(ll)€lqr£lq8rmdS'k 

J s 
+ J Cij k iGjnMK){bi+€ lqt n q (xl-x;)}dSi n 

^(p €pmkCij k iGj n j(R){bi-\-ei qr Q, Q (x' r — x° r )}dL' p 

+ I ^p>nkCi jk iGj» i i(R)Cl p dS' l + p* n (r). ( 534) 

•'iS 

Here the second equality follows from Stokes' theorem (A4), and the third from (2.7), (B5), and (5.12). 
Now if we call /3,„„ = u T n m — fi* H the "elastic distortion," it is not a state quantity because the 
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second integral in (5.34) cannot be written as a line integral. The proof of this consists merely in 
showing that the integrand in this term is not divergence-free; using (2.7) we find 

\€pmkCijklGjn,i{r)\ ,(= €pmkC ijklGjn ,n(r) 
= — €pmn&(r) 

# 0. (5.35) 

If we compare the second integral of (5.34) with (3.30) we notice a great deal of similarity. As is 
well-known the jump of u T n in (3.30) across the surface S is given by \u^\ = b n according to (3.18). 
From this we deduce that the second integral in (5.34) leads to a jump of the total distortion f} T mn 
across the surface S of [jS^J = e P mn£l P - This is a jump in rotation of [co£] = ll2e P mn[fi T mn \ = Qp, 
in agreement with (5.6). In the older approach to dislocation theory where the surface S was ignored, 
the displacement of a discrete dislocation line was regarded as a multiple-valued function with a 
period of the Burgers vector. Equation (5.34) shows that from this point of view the rotation (or 
distortion) of a discrete disclination line is a multiple-valued function with a period of the Frank 
vector. 

The elastic strain is from (4.14), (5.14), and (5.34) 

e mn (r) = <b € P mkCijkiGj n ,i(R){bt+ ei qr £l q (x' r — x°)}dL p{mn) + I €pmkCij k iG jn j(R){l p dSl (rnn) . (5.36) 

The integrand of the second integral now is divergence-free and so it can be written as a line 
integral 

— I eq8lImnpq,8(R)Q>pdS' l =— <b I mn pq(R)£lpdL q 

by (4.28) and Stokes' theorem. Hence 

e m n(r) = <b ep mk C m G jn j(R){bi^ei qr n q (x f r -xl)}dLl imn) -i> I mnpq (R)Q,pdL' p . (5.37) 

This result could of course also be obtained directly from (4.29), (5.19-20) and (Bll). This is the 
relation we sought. We see that the elastic strain can be written as a line integral along the discrete 
defect line. Therefore it is a state quantity. 

Equation (5.37) also allows us to estimate the asymptotic behavior of the strain at large dis- 
tances from a straight disclination line. The second integral will give the dominant term. The 
incompatibility source tensor Imnpq(r) varies as r _1 . Due to the integration then e mn (r) will vary 
as In r as r— > oo. Since the displacement is an integral of the strain, it will vary as r In r as r—> °°. 
These results are also listed in table 3. 

We wish next to derive the elastic bend-twist. First we find the derivative of the total distortion 
(5.34) 

u Z,ms( r ) = ( f ^pmkCi jk iG jn ,is(R){bi-\-€i q ril q (x , r — x r )}dL , p -\- I € P mkCwaGj n ,u(R)ClpdS' i + j8* n> ,(r). 

(5.38) 
The second term becomes by Stokes' theorem 

i> €pmk€qslCijklGjn t i(R)flpdL q -\- J €pmkCijklGj n ji(R)Cl p dS f s . 
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The second term in this expression becomes by (2.7), (B5), and (5.13) 

— €pmkdkn&(R)ttpdS' a = ~ €pmn8 a (S)(l,, = €/>//,,,</>*, (r) . 
Hence (5.38) becomes 

uT n,ms^ = J CpmkCijklGjn, is(R) {bl+ €l qr £l q (X r — X° r )}dL' p 

+ 4> €pmk£qslCijklGjn , i ( R ) £lpdL ' q (5.39) 



+/Cjw+*«*;(r). 



We now find from (4.15) and (5.15) 

Kgt(r) = 1/2 <J) etmnepmkCijkiG jn ,i 8 (R){bi + eiqr£l q (x' r — x*)}dL p 

+ 1/2 (1) €tmn€pmk£qslCijklGjn, i(R)ClpdL' . 



(5.40) 



This could also have been obtained directly by substituting into (4.34) from (5.19-20). It is also seen 
to be a state quantity. 

This section has extended the results for a discrete dislocation line of section 3.2.2 to a more 
general discrete defect line. The central result there, equation (3.32), has been replaced by (5.37) and 
(5.40). Again we see the elastic distortion of dislocation theory is generalized into the basic elastic 
fields of defect theory. 

6. Continuous Distribution of Infinitesimal Defect Loops 

In section 4.1 the defects were defined by (4.1-2) and the characteristic vectors by (4.5-6) in 
terms of the given basic plastic fields, strain and bend-twist. Eliminating the basic plastic fields from 
these definitions led to (4.7-8), relations between the characteristic vectors and the defect densities. 
These equations could alternatively be used to define the defect densities in terms of the character- 
istic vectors, if they are prescribed in a suitable manner. It is convenient to put these relations into 
differential form. First we define the disclination density. When the defects are continuously distrib- 
uted the disclination density tensor is defined locally by 

Opq-js;- (6.1) 

For a distribution of discrete defect lines this represents the average disclination density where ACl q 
is the qth component of the resulting Frank vectors of all the disclinations which pierce through a 
surface element AS P oriented normal to the x p direction at the given point. In view of (4.7) it is neces- 
sary to modify the definition (3.33) for the dislocation density a p i. For an inhomogeneous continuous 
distribution of defects we define it locally by the relation 

_ &Bi 

OLpl — eiqrOpqXr = To~ * (6.2) 

For discrete defect lines this relation will give us the average dislocation density, where ABi is the 
/th component of the resulting total Burgers vectors of all the defects (dislocations and disclinations) 
which pierce through a surface element AS P normal to the x p direction at the given point. 
Reference to (5.31) suggests that a more convenient way to write (6.2) is 



AW 
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apl = ^s7' (6 - 3) 



where we have defined 

A [Ui] = A#, + eiqrkClqXr, (6.4) 

as suggested by (5.1). For discrete defect lines A [ui] represents the /th component of the resulting 
displacement jumps of all the defect lines which pierce through a surface element AS P oriented 
normal to the x p direction at the given point. Relation (6.3) is what would be expected if the dis- 
location density is determined by the local plastic deformation, regardless of whether it is due to 
discrete dislocations or disclinations. For generality relation (3.33) could also have have been put 
in the form (6.3) by reference to (3.18). 

We next consider a continuous distribution of infinitesimal defect loops differently oriented in 
space. We introduce the disclination loop density tensor following Kroupa's [5] line of reasoning. 
We can define it as follows: (# q represents the flux of disclination (or Frank vector) in the ^direction 
that encloses a unit vector in the Xu direction. When the loops are continuously and inhomogene- 
ously distributed the density tensor t^q is a function of the position and is defined locally by 

_ Afl g 

^ = aZ7 (6 - 5) 

For a distribution of discrete loops this represents the average disclination loop density where now 
Aflg is the q\h component of the resulting Frank vectors of all the loops which are pierced by the 
line element AL/> oriented in the Xk direction at the given point. It is now also necessary to modify 
the definition (3-34) for the dislocation loop density. For a continuous distribution of defect loops 
it is defined locally by 

y»—fiZ~ (6-6) 

For discrete loops this represents the average dislocation loop density, where now A[w/] is the /th 
component of the resulting displacement jumps of all the loops which are pierced by the line 
element AL A oriented in the x k direction at the given point. By (6.4) we also have the alternative 
definition 

Aft, All, 

yki =u: k + eiqr ^: k Xr ' (6 - 7) 

where for discrete loops Aft is the /th component of the resulting total Burgers vectors of all the 
defect loops which are pierced by the line element AL/c oriented in the Xk direction at the given 
point. 

To derive the relations between the defect loop densities and the basic plastic fields, we first 
combine (6.5) and (6.7) into 

Aft 

Jkl — tlqrCkqXr = TjT ' (6.8) 

The relations (6.8) and (6.5) are easily converted to integral form 

Bi = j> (yki — eiqrbcqX r )dL k9 (6.9) 

CkqdLk. (6.10) 



VL q = j> 



Now we compare these relations with the definitions (4.5-6). The integrands can only differ by a 
gradient with respect to x k > Therefore we can set 
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kq 



-&«+<*, (6.11) 



e kl~~ € lQr K kq x r=—yi x i+€i Qr £ > k<pCr+ (tf — €lqr<0 P JCr) Jc, (6.12) 

where w 1 ' and u l * are arbitrary vector fields, subject only to the condition that e p is symmetric. This 
last condition will provide a relation between them, as we shall show. If we substitute (6.11) into 
(6.12) we obtain 

««*-?«-*«*£+ "(!*• ( 6 - 13 ) 

Now the symmetry condition on e p t gives the relation between o) 1 ' and u p 

<tf q = -\l2€ k i q {y k i-u p lk ). (6.14) 

From this we get for (6.13) and (6.11) 

<i = -y<*0+<*), (6.15) 

K mq = ~ {»«~ Il2€ k i q y k i, m +ll2€ k iqii p fcm . (6.16) 

These are the relations that identify the basic plastic fields for a continuous distribution of defect 
loops. With them we can find all the relations derived in section 4 in terms of a loop distribution. 
For example, from (4.1-2) we find 

Otpl= € P mk{ykl, m + €klq£>rnq) , (6.17) 

&pq = €pmk£>kq, m« 

(6.18) 

This is the fundamental relationship between the defect loop densities and the corresponding 
defect densities. We see from (6.17) that the disclination loop density (, mq also contributes to the 
dislocation density, in addition to the dislocation loop density y^b 

Since a plastic displacement such as u p in (6.15-16) does not contribute to the elastic fields 
(c.f. section 2.6), we can set u p =0 without loss of generality, and so we can use 

* = -ycw>, (6.19) 

K mq = -£>mq-ll2e k iqykl, m (6.20) 

for the purpose of calculating the fields of a continuous distribution of infinitesimal defect loops. 
Furthermore, we note from (2.15) and (6.19) that only the dislocation loop density y^i will contribute 
to the total displacement. Hence the elastic strain and stress are unaffected by the disclination 
loop density, £ m(/ , as we already hinted at in section 5.3. We also note that all the above results 
reduce to those of section 3.3 when the disclination loop density vanishes, J m? = 0. 

Mura [4] generalized his "plastic distortion" and "plastic rotation," which he had defined for 
a discrete loop as in section 5.2, to a continuous distribution. We shall now interpret his approach. 
If we compare (5.14-15) with (6.19-20) we can make the following identification: 

^mq + g £ *W, m = ~ £mq ~ CkhfYkl, m- (6.22) 
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These equations can be solved for Mura's plastic quantities as follows: 

ffii=-yki-€ki**^ (6.23) 

<£m 9 =-£m 9 +< m , (6.24) 

where w* is an arbitrary vector field. These relations identify Mura's quantities for a continuous 
distribution of defect loops. The basic plastic fields are obtained in terms of Mura's quantities by 
substituting in (6.19-20) 

*=/B£* (6-25) 

< Q = Kq + ^klqPtt, rn. (6-26) 

These relations are identical in form with (5.14-15), but are now also valid for a continuous distri- 
bution. We find the defect densities by substituting in (4.1-2) 

otpi = — e pmk (j8* f , m + e k iq(t)Zq) , (6.27) 

6pq = ~ €pmk<i>kq, m, (6.28) 

which are identical to relations in section 5.2. The interesting vector o>g does not affect the basic 
plastic fields and the defect densities, and hence it does not affect the elastic fields either. In (6.23) 
it contributes the antisymmetric term — €ki q (Oq to /3^. Mura [30] has called the dislocations resulting 
from an antisymmetric plastic distortion an impotent distribution of dislocations, because, as can 
be deduced from (2.15), such a distortion does not contribute to the total displacement, and hence 
gives no elastic fields. In general an antisymmetric plastic distortion will give a finite dislocation 
density, c.f. (3.4). However, the term a)* m in (6.24) is exactly right to annihilate both defect den- 
sities (6.27-28) due to (Oq. Hence, we can set a)q = without loss of generality for the purpose of 
calculating the elastic fields due to a given distribution of /3*, and</>* g , or 

fi?i=-y«. (6-29) 

This shows that Mura's "plastic distortion" and "plastic rotation," introduced in section 5 for a 
discrete defect line, can be interpreted as the dislocation and disclination loop densities, except 
for a minus sign. This then resolves a difference we had with Mura. Relations (6.25-26) are the basic 
relations that connect Mura's approach with ours. For example the characteristic vectors are 
found from (4.5-6) to be 

Bi=-&(Pi-6 k udtjCr)dL k , (6.31) 

Cl q =-j^* q dL k9 (6.32) 

which correspond to (6.9-10). 

We can now also give an interpretation to (5.12-13) in terms of infinitesimal loops. To construct 
a discrete defect line L, we first distribute a constant density of infinitesimal disclination loops of 
strength [\ Q over any surface S whose boundary is L. This distribution only gives a rotation across 
S, but no stress, and furthermore no unique axis is defined. Rather, each infinitesimal loop has its 
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own axis, so that the resultant axis is continuously distributed over S. Now we add a distribution 
of infinitesimal dislocation loops to S, consisting of two parts. One, containing the dislocation 
Burgers vector, is a constant distribution of strength bi, which gives the discrete dislocation line. 
The other, containing the Frank vector Q, Q , has just the right strength to bring the axis of each 
infinitesimal disclination loop to the point x° r . We see that this is done by increasing the dislocation 
loop strength proportional to the distance from the chosen axis. The combination of this linear 
infinitesimal dislocation loop density with the constant infinitesimal disclination loop density men- 
tioned above gives the discrete disclination line. It is the infinitesimal dislocation loop distribution 
that gives rise to the elastic strain of the discrete disclination line. 

So we conclude that for a finite defect loop the defect loop densities are given by 

y k i(r) = d k (S){b l + 6i Qr n Q (x r -xO r )}, (6.33) 

Ck q (r) = 8 k (S)Sl q . (6.34) 

The vector 8k(S) is the Dirac delta function on the surface S and it is always normal to S. In a 
plane we have a twist or a wedge disclination loop according to whether the Frank vector is normal 
or parallel to S, respectively. Therefore (6.34) shows that the diagonal and off-diagonal components 
of £a-<7 represent the twist and wedge components of the disclination loop density, respectively (see 
table 1). Equations (6.33-34) show how to make the transition from a continuous distribution of 
infinitesimal defect loops to a finite defect loop. By (B15) and (5.17) we see that they satisfy (6.9-10), 
which also remain valid for many finite loops, and therefore can be used for the transition from finite 
loops to a continuous distribution: For many finite defect loops the average dislocation and disclina- 
tion loop densities jki and £a- 9 represent the sum ol the x\ component of the displacement jump 
vectors and the x q component of the Frank vectors, respectively, of all the loops whose surfaces 
are intersected by a unit vector in the xa direction. 

The infinitesimal defect loop density tensors ^q and yki were defined in this section independ- 
ent of the definitions of the defects densities in section 4.1 and the discrete defects in section 5.1. 
What is the correlation? 

Let us first examine the relation between continuous distributions of defects and loop den- 
sities. For dislocations only there is a straightforward relation between them given by (3.37). When 
disclinations are introduced there is a similar straightforward relation between the disclination 
loop density and the disclination density tensor, given by (6.18). However, the dislocation density 
now contains a contribution from the disclination loops, equation (6.17). This means that the two 
definitions of the "dislocation" are not identical, and that a certain amount of mixing occurs in 
going from one formulation to the other. 

Now let us examine the relation between a discrete defect line and a continuous distribution 
of loops. For dislocations only there is again a straightforward relation between the dislocation loop 
density and the Burgers vector, given by (3.39), or between the Burgers vector and the dislocation 
loop density, given by (3.35). When disclinations are introduced there is a similar straightforward 
relation between the disclination loop density and the Frank vector, given by (6.34), or between the 
Frank vector and the disclination loop density, given by (6.10). However, the dislocation loop den- 
sity for a discrete defect line now contains a contribution from the Frank vector, equation (6.33), 
or the general Burgers vector contains a contribution from the disclination loop density, equation 
(6.9). Therefore the definitions of the dislocation are not identical in both approaches, and as we see 
from the relations quoted, a certain amount of mixing occurs in going from one formulation to the 
other. 

There are, therefore, at least three independent ways to define the dislocation content of de- 
fects: in terms of a continuous distribution, a discrete line, or a continuous distribution of infini- 
tesimal loops. For dislocations only, these definitions are equivalent but with disclinations they are 
essentially different. 
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7. The Discrete Dipole Line 

7.1. Basic Relations: The Biaxial Dipole 

Extending Kroupa's [18] definition of the dislocation dipole, we define the discrete dipole 
line as a close pair of discrete defect lines with opposite characteristic vectors. We call it the 
biaxial dipole because this defect would have two rotation axes. We wish to present in this section 
some of the relations analogous to those for the dislocation dipole in section 3.4.2. From (3.42), 
we find the defect densities of the biaxial dipole conjugate to the basic dislocation and disclination 
densities and 6 pq to be 

o$i(r)=-€ma p i,m(r) 9 (7.1) 

0£ 9 (r) = -£n0p*,m(r). (7.2) 

By (5.19-20) these relations become for the discrete dipole line 

<(r)=-8 P , m (L){6/ + €/ gr a <7 (x r -<)}^-Op(L)€^a^ r , (7.3) 

0&(r)=-8 P , m (L)ft«6n. (7.4) 

The meaning of the displacement g m is as follows: We let £m~" * 0, 6|— * <», and £l q -> °° in such a 
way that b^ m and Q,qg m remain constant. 

The displacement of the dipole line is from (3.44) 

"£(r) = -6< m (r) (7.5) 

= ~j> e P mkCijkiGj„, ,(R) {b, + ei qr Clq(x' r - x° r ) }^ m dL' p 

I GpmkCijklGjn, j ( R ) fl,,f m<£SJ 

+ 8 m (S){b l + e„ q ril q (xr-x°)}{ m , (7.6) 



■{. 



where we have used (5.34) and (5.12). In this expression the third line can be made to vanish by 
choosing S near the line to lie along £, or g m dS' m = along L. Nevertheless, due to the surface 
integral, the displacement is not a state quantity in contrast to the case of the dislocation dipole. 

From (7.1-6) all other relevant quantities for a discrete biaxial dipole line can be derived if 
desired. 

Equation (7.6) allows us to estimate the asymptotic behavior of the displacement at large 
distances from a straight biaxial dipole line. It is easily deduced that it will vary as In r as r^ o°. 
Hence the strain will vary as r _1 as r— > oo. These results are also listed in table 3. 

7.2. Influence of the Axis: The Uniaxial Dipole 

The dipole of section 7.1 is obtained from the basic discrete defect line of section 5 by trans-, 
lating it through the infinitesimal distance f, including the axis of rotation, and then subtracting 
the basic defect at its original position. 

In this section we first wish to isolate the influence of the motion of the axis, when the position 
of the defect line is held fixed. We therefore consider the following defect, again composed of two 
parts: the first is obtained from the basic defect by translating its axis only through an infinitesimal 
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distance £, and the second part is the negative of the basic defect at its original position. The 
defect densities of this defect are given by 

af^in.da^ldxl, (7.7) 

By (5.19-20) these relations become 

c$ l (r) = -8 P (L)ei v n g gr, (7.9) 

<&(r)=0. (7.10) 

The displacement is 

^(r) = f m a^(r)/d^ (7.11) 

= ~" Jr CijklGjn, i(R)eiq r Q q £;rdS' k 

= — I epmkCfjkiGjn, i(R)n p g m dS' r (7.12) 

where the second equality follows from (5.33), and the third by rearranging the indices using the 
symmetry condition (2.3). 

n„ ~~, : — tn a\ 1 <n io\ ,.,:*u to oo\ „~J /o om ,..„ *u~* *l . _♦ *L~ J:~l *:~~ 

wn tuuipainig, \t.sj aim \^.±^) vvilii \v«^*-V anu \*j.ij\j), v»l- oi^<o tuat luc-j it/^icocut iin., uiaju^anun 

density and displacement of a discrete dislocation line with Burgers vector 

bi=-e, qr ft^r. (7.13) 

This was to be expected, since the motion of the axis through the distance % r has the effect of 
translating the two sides of S by the distance bi given in (7.13), as we can see from (5.5). This con- 
clusion complements the statement in section 5.2 that the axis of a discrete defect line can be 
moved by adding a dislocation to the line, c.f. (5.23). 

Next we want to examine the effect of holding the axis fixed and moving the defect line only. 
This type of the defect will be called a uniaxial dipole because it has only one axis. It could alterna- 
tively have been used as the definition of the discrete dipole line. It is composed of the following 
two parts: the first is obtained from the basic discrete defect line by translating it through an infin- 
itesimal distance £ keeping its axis fixed, and the second part is the negative of the basic defect 
at its original position. The resulting fields are simply the difference between those in section 7.1 
and the above. So the dislocation density and displacement of the uniaxial dipole are 

a L pl {r) = a° l {r)-a%{r) (7.14) 

= -d p , m (L){b,+ e, q rn q (xr-x<!.)}£ m , (7.15) 

<(r) = <(r)-<(r) (7.16) 

= -i e pmk Cij kI Gjn,i(R) {bi + ei Qr n q (x' r -x° r )}t m dL' p , (7.17) 

whereas the disclination density is the same for either type of dipole, 

6' m =6%. (7.18) 
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We note that the displacement of the uniaxial dipole is a state quantity. 

A special case of (7.17) is of particular interest, namely, the one corresponding to the wedge 
disclination coinciding with its axis as the basic defect. In this case 6/ = 0. For a wedge disclination 
the Frank vector is parallel to the disclination line, ft||dL'. The point r' lies on the line and the 
point r° on the axis. Since the line coincides with the axis, the difference, is also parallel to the line, 
(r'-r°)||dL\ So ft||(r'-r°), and hence ei qr (l q (x f r -x° r ) =0. We conclude that (7.17) vanishes 
for this case. Therefore u% = u*. This means that the conjugate dipole, corresponding to this wedge 
disclination with Frank vector £l q , is the dislocation with Burgers vector bi given by (7.13). It is an 
edge dislocation because b _l_ il||dL\ Eshelby [21] used this approach to give a simple derivation 
of the elastic field of an edge dislocation, when the field of the wedge disclination is known. 

Equation (7.17) shows that the asymptotic behavior of the displacement at large distances from 
a straight uniaxial dipole is the same as for a biaxial dipole. Table 3 shows there is a gap in the 
asymptotic behavior between a disclination dipole and loop, i.e., there is no disclination type defect 
with the r~ l behavior for the displacement. 

8. Application to Dislocations 

The purpose of this section is to check the internal consistency of the results for dislocations 
and the more general defects. It contains no new material. We want to show that the results of 
section 3 fall out of section 4-6, when no disclinations are present, i.e., when the plastic deformation 
of the body is completely described by the plastic distortion jSj^. 

8.1. Continuous Distribution of Dislocations 

8.1.1. Geometry 

The plastic strain and rotation are now given by (c.f. 3.2-3) 

*Zi = P(kih (8.1) 

<=l/2e^/%, (8.2) 

so that 

e& + €*^=/3&, (8.3) 

in agreement with (3.1). The plastic bend-twist is (c.f. 3.12) 



C kq W 9,A- 



(8.4) 



Hence we find from (4.1) 

OLpl = — Cpmk(e% lfm + eklqto^m) =~ €pmkPkl,m> (8.5) 

in agreement with (3.4), and from (4.2) 

Opq = — tpmktoq, km = ®> (8.6) 

as expected. The Burgers vector is found from (4.5) 

Bi = —$> (e^ — eiqrio^jpc^dLk 

(8.7) 



= -<f (3 p kl dL k . 
A 
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Here the first equality follows from (8.4), the second from a partial integration where the integrated 
part vanishes around the closed curve X, and the third from (8.3). The result agrees with (3.7). 
Furthermore we also see that (8.6) in (4.7) agrees with (3.8). The Frank vector is found from (4.6) 
and (8.4): 



ilq = — <j> (Oq,kdLk = 0, 



(8.8) 



as expected. This also agrees with (8.6) in (4.8). 

In a similar way it is easily shown that (4.17) and (4.20) reduce to (3.10) and (3.11), whereas 
(4.18) and (4.21) vanish. 

8.1.2. Statics 

The Kf p term in (4.23) becomes by (8.4) 

— I e P mkC i j kl C7jn,i{K)o) I l )ll {r')dV'=- I z pmk C y**Gj n ,*/(R)G>£(r')dF' 

= €pinnO)p(r). 

Here the first equality follows by a partial integration, and the second by (2.7). Hence we have for 
the total distortion (4.23). 

u Ti,m( r ) = I ^pmkCi jk iG jn j(R)a v i(r f )dV'^€^ n (r)^€t >mi iOJ^(r), (&Q) 

in agreement with (3.16) and (3.1). From (4.14) we find the elastic strain 

€mn(r)= I € P mkCijklGjn t i(R)<Xpl( r ')dV' (mn). (8.10) 

This relation can also be obtained from (4.29) with (8.6). It is in agreement with (3.13) and (3.17). 
We next find the bend-twist. From (8.9) and (8.4) we find 

u l,msM = I ^pmkCijkiGjnM^oipiir^dV' ^ e^ nns (r) -^ ep m nK^ p (r) . (8.11) 

Hence we find from (4.15) 

K st (r) = 1/2 J €tmn€pmkC m £ in AB)a p i(v')dV' . (8.12) 

This relation can also be obtained from (4.34) with (8.6). We see it is in agreement with (3.14) and 
(3.17). 
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8.2. The Discrete Dislocation Line 

8.2.1. Geometry 

When a discrete defect line contains no disclination, Cl Q = 0, as we saw from (8.8). For this 
case (5.5) reduces to (3.18) and (5.6) becomes 

[a)*(r)] = 0, (8.13) 

while (5.7) reduces to (3.19) and (5.8) reduces to (3.20). 
Next we find that (5.12-13) become 

/B&(r)=-MS)6i, (8.14) 

<(r)=0. (8.15) 

This shows from (3.21) that 

dEi=/B&. (8-16) 

as expected, since both sides represent the dislocation loop density by (6.29) and (3.38). Equation 
(8.15) shows that the disclination loop density vanishes. Now (5.14-15) become from (8.14-15) 

e? l (r)=-ll2[8 k (S)bi+8i(S)b k ], (8.17) 

K mq( r ) = ~ V2Cklq8k,m(S)bi, 

(8.18) 

(K^(r)=-H2e lqr drAS)b l ). 

We see that (8.17) is in agreement with (3.2) and (3.21) and that (8.18) agrees with (3.12) and (3.21). 
From (5.17) we find that 

B, = b h (8.19) 

showing that the total Burgers vector reduces to the dislocation Burgers vector. From (4.1) and 
(8.17-18) we find 

<x p i(r)= € pm k8k,m(S)bi= 8 p (L)bi. (8.20) 

This result can also be obtained directly from (5.19) and is in agreement with (3.23). Next we find 
from (4.2) and (8.18) or directly from (5.20) 

Mr) = 1/2 e pmk eiqr 8 r , km (S) bi = 0, (8.21) 

as expected, in agreement with (8.6). 

8.2.2. Statics 

Again for Cl q = 0, we see that (5.33) reduces to (3.30), and (5.34) reduces to (3.31) with (8.16). 
The elastic strain is 



e mn (r) = (p epmkCijkiGjn, i(R)bidL' 

J j r (mn)> 



(8.22) 



This result can be found in two ways: first from (5.37), or second from (4.29) with (8.20-21). It is in 
agreement with (3.13) and (3.32). 

Next we find the elastic bend-twist from (5.40) 

K. v ,(r) = 1/2 d) € t mn€p,nkC^kiGjn, i s (R)bidL' p . (8.23) 

This could also have been found from (4.34) with (8.20-21), and is in agreement with (3.14) and 
(3.32). 

So we have shown in this section that the more general defect theory, including disclinations, 
completely reduces to the well-known dislocation theory in the special case that the disclinations 
vanish. 

9. The "Dislocation Model" of a Discrete Defect Line 

Li and Gilman [20] considered the finite disclination loop as a continuous distribution of dis- 
locations, and called this the "dislocation model" of the disclination. Mura [22] used the same 
concept, which he ascribed to Eshelby, discussing also the case where he replaced a wedge dis- 
clination by a semi-infinite edge dislocation wall. We wish to make clear the distinction between the 
two concepts. 

We start with the observation, made in section 6, that only the dislocation loop density, or 
Mura's "plastic distortion" /3^, contributes to the elastic strain. For a discrete defect line, this 
quantity is given by (5.12). So the elastic strain obtained from this expression does not depend on 
what we choose for the disclination loop density, or Mura's "plastic rotation" c/)*^ . If we choose it 
to vanish 

Ki = °> (9.1) 

then we have a distribution of dislocation loops over the surface S that gives exactly the same elastic 
strain as the discrete defect line of section 5. We call the corresponding dislocation distribution the 
dislocation model of the defect line. We note that this dislocation model is clearly a different defect 
from the discrete defect line it corresponds to. Specifically, we obtain the dislocation model by 
setting 

££=#„ (9.2) 

where 

/3\i(r) - - 8 k (S){bi + e lqr n Q (x r -X« )} , (9.3) 

and where bi and £l Q are constants. We can then use the methods of section 3.1 to find any other 
desired relations. Furthermore (5.14-16) become by (9.1) 

< = 1/26*1^,,!, (9.5) 

*«.m+€*iflici^=j8?| -m . (9.6) 

With these relations we can alternatively use the methods of section 4 to find any other desired 
results 
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9.1. Geometry 

We find the dislocation density from (3.4) with (9.2) or (4.1) with (9.6) 

Ctpl{r)=-€pmkfikl i m 

= €pmk&k,m(S) {b I + €iqr£lq(x r ~ X r ) } + € pm k$k (S)6( qr £l (firm 

= 8 p (L){bi+ci^ q (xr-x«)} + 8i(S)a p -8pfi k (S)a k . (9.7) 

Here the second equality follows from (9.3), and the third from Stokes' theorem (B26) and the 
expansion (A3). We see now that the dislocation model is a dislocation distribution consisting of 
the following two parts: First a dislocation line along L, which is the same as for the corresponding 
discrete defect line, (5.19), and second a constant dislocation distribution over the surface S. In 
other words it is a dislocation wall at the surface S, which terminates on a dislocation line at its 
boundary L. The disclination density vanishes from (4.2) and (9.5) 

PQ =-e pm kll2e lQrr 3Xkrn = ^^ (9.8) 

as expected. 

We see that the continuity equation (3.6) or (4.3) is satisfied 

<*pi, P (r) = *p.p(L) {bi+ € t<p Sl q (x r -x° r ) } + 8 p (L)€i qp Sl q + Bi, p (S)a pr 8*,,(S)n* 

= 0, (9.9) 

where we have used (B27-28). 

We next find the total Burgers vector from (4.5) with (9.4-5) 



5/ = 



j> x <fi5 a >rPlri\.*Xr)dL k 

pkidL k 

= j> 6 k (S){b l ^€i qT n q (Xr-X r )}dL k 

= bi+€ lqt n q (x s r x -x° r ). (9.10) 

Here the second equality follows by a partial integration, the third from (9.3), and the fourth from 
(B15-16), where xf x is the point where X crosses S (fig. 1). 

9.2. Statics 

If we substitute (9.3-4) into (2.15) we find for the displacement 

u T n (r)= f C ijk iG jn>i (R){bi + 6 lq rn q (x r -x r )}dS r , (9.11) 

which is identical to (5.33). So we see that the displacement of the dislocation model is exactly the 
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same as the displacement of the corresponding discrete defect line. Hence (5.34) follows in exactly 
the same way 

Uljr) = * epmkC V kfijn,l<JR){bi + e lq M„{x' r ~ X° r )}dL' v + J epmkCvkGjn^tlpdS', + £*>)•• 

(9.12) 
So we find the elastic distortion from (3.9) and (9.2) 

Pmn(r) = <j) epmhCi jk iG Jn ,i(R){bi + €l qi £l q (x' r — x) }dL'+ I €pmkC klmn G jn , i(R)ft pdS\ (9.13) 

This relation could also be obtained from (3.17) and (9.7). It follows from the discussion in section 
3.1.2 that this elastic distortion is a state quantity for the dislocation model, but according to the 
discussion of section 5.3 it is not a state quantity for the discrete defect line. This may be the physical 
significance of Mura's "elastic distortion" for a discrete disclination. 



emn(r)=6€pmkCijkfijnA^){bl^€lqrCl q (xl-xl)}dLp (mn) 

+ I e P mkC ij k iG j?l ,i(R)fi p dS f , (mn) . 



(9.14) 



This relation also follows from (4.14) with (9.4) and (9.12), or (4.29) with (9.7-8). We see that this 
expression is identical to (5.36). So the elastic strain and hence stress of the dislocation model is 
identical to that of the corresponding defect line, as we stated before. 
We next find the bend-twist. From (9.12) we have 

/• 
a n,m«( r ) = ( P epmkCi jk iGjn,is(R){bi-\-e lqt n(x' r — x r )}dL p 

(9.15) 
+ J € pmk CijkiQjn t is(R)il p dS , l +Pl nt8 (r) 9 

and from (4.15) and (9.5) 

K st (r) = 1/2(1) €t m n€pmkCij kl G jn j s (R){bi + ei q rU q (x' r — X r )}dL p 



1/2 J €tmn€pmkGjn,is(R)npdS' l . 



(9.16) 



This relation could also have been obtained from (3.14) with (9.13), or from (4.34) with (9.7-8). 

We see from this section that there is a great similarity between the dislocation model and its 
corresponding discrete defect line. Therefore it is important to distinguish carefully between them. 
For example, the defect densities for the discrete line are given by (5.19-20), while for the dis- 
location model they are given by (9.7-8). So we see that in the transition the disclination density 
in (5.20) has been traded for the constant surface dislocation density in (9.7). As we saw both 
defects give the same elastic strain; explicit expressions for it can be obtained by substituting the 
above densities into (4.29). For the discrete defect line we find (5.37) and for the dislocation model 
(9.14) for the elastic strain. We see that the second term in (5.37) represents the contribution from 
the disclination density; it equals the second term in (9.14) which represents the contribution from 
the constant surface dislocation density. We shall illustrate some special cases of the dislocation 
model in future publications [25, 26]. 

Li and Gilman have also calculated the force on a discrete disclination line by assuming that 
it equals the force on the corresponding dislocation model. Since the present paper does not deal 
with forces a formal proof of this interesting result will be published elsewhere [23]. 
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10. The Compensated Disclination Line 

We now want to investigate a problem that complements the one of section 9. Consider a 
vanishing dislocation loop density and a disclination loop density given by (5.13) 

0* = O, (10.1) 

*£(r)=- MS)O t , (10.2) 

where Q q is a constant. This represents just a constant distribution of disclination loops on the 
surface S. Since there is no dislocation loop distribution, we shall find that there is no elastic 
strain for this defect. We find the plastic quantities from (5.14-16) 



oi y -. 

'kl 



< q =K- 



(10.3) 
(10.4) 



We now use the methods of section 4 to find any further desired results. 

10.1. Geometry 

The dislocation and disclination densities are from (4.1—2) 

OCpi(r) =e P mk€klq$m{S)£l q 

= 8 pi 8 A .(S)ft fr -8KS)n p , (10.5) 

0p Q (r)=6pmkdk,m(S)n q 

= 8p(I)n ff (10.6) 

where we have used Stokes' theorem (B26). Therefore the defect in this case is a dislocation wall 
at the surface S, which terminates on a disclination at its boundary L. It consists of the same dis- 
crete disclination line as treated in section 5, (5.20), and a dislocation distribution on S, which is 
just right to make the elastic strain vanish. Therefore we have called it the compensated disclination 
line. Note that the sum of (10.5) and (9.7) gives (5.20), and that (10.6) is the same as (5.19). This 
was of course to be expected, because sections 9 and 10 represent a decomposition of the problem 
of section 5. We find that the above results satisfy the continuity equations (4.3-4) 

a pl , p (r)^€i PQ 6 I)Q (r) = d k ,i(S)n k -d l , p (S)n p ^€i pq d p (L)CL Q =0, (10.7) 

pq , P (r) = 8 p , p (L)n q =0, (10.8) 

using (B27-28). 

We next find the total Burgers vector from (4.5) and (10.2^4) 



Bl = — <j> €l Q rdk(S)Ct q X r dL k 



= -€, qr n Q x s r \ (io.9) 
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by (B16), where xf x is the point of intersection between the curve A and the surface S (fig. 1). We 
note that the sum of (10.9) and (9.10) is the same as (5.17). We find the Frank vector from (4.6), 
(10.4), and (10.2) 



-j> K* q dL k = j> 8 fc (S)n^ fc = n g , (10.10) 

using (B15). This identifies the constant £l q as the Frank vector. The same results can be obtained 
from (4.7-8) and (10.5-6). 

10.2. Statics 

If we substitute (10.3) into (2.15) we find the displacement 

<=0. (10.11) 

Hence from (4.14) and (10.3) 

e fc *=0. (10.12) 

So we see that there is no elastic strain and hence stress, as we discussed before. The elastic bend- 
twist is from (4.15), (10.11), (10.4), and (10.2) 

For consistency we show that this relation can also be derived another way. If we substitute (10.5-6) 
into (4.34) we find 

Kst(r) = — 1/2 letmnejnnkCijklGjnjsi^CLpdS) + 1/2 (I) €tmn€pmk€q8lCiJk£jn, i(R)£lpdLq 9 (10.14) 

which by Stokes' theorem (A2), (2.7), and (B5) becomes 

fCtf(r)=— 1/2 I €tmn€pmkCijklGjn,il(R)0,pdSg 
= 1/2 I €fmn€pmA-8fcn8(R)ftp(£S' g 

= 8 8 (S)a t) (10.15) 

in agreement with (10.13). Note that the sum of (10.14) and (9.16) gives (5.40), as expected. 

Summarizing, we see that when the compensated disclination line is added to the dislocation 
model we obtain the discrete defect line. 

11. Summary 

We started this paper with the general solution of the plastic strain problem which is just an 
extension of Eshelby's transformation problem and essentially equivalent to Mura's plastic dis- 
tortion problem. It formed the basis of all static defect fields. We then reviewed dislocation theory, 
including the continuous distribution, the discrete line, Kroupa's continuous distribution of in- 
finitesimal loops, and the dipole. This introductory material formed the point of departure for the 
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general theory of defects, i.e., disclinations and dislocation combined. On the other hand, it formed 
a basis of comparison, because the general defect theory reduces to it when the disclinations 
vanish. 

We defined the continuous distribution of defects. This definition was motivated by a violation 
of the compatibility equations. We also defined the Frank vector, which is the characteristic 
rotation vector of the disclinations, analogous to the Burgers vector for dislocations. We derived 
closed integral expressions for the basic elastic fields, the elastic strain and bend-twist, in terms 
of the defect densities, showing that they are state quantities. These integrals contain kernels 
with Green's tensor and the incompatibility source tensor, a type of Green's tensor introduced 
by Simmons and Bullough. 

The definition of the discrete defect line was motivated by Weingarten's theorem. A new 
quantity introduced here is the axis of rotation, which did not exist for a continuous distribution. 
After finding the basic plastic fields, the plastic strain and bend-twist, all the results for the con- 
tinuous distribution can be specialized to the discrete case. We found that the calculations were 
simplified by using two new quantities, Mura's "plastic distortion" and "plastic rotation," which 
we later identified as the dislocation and disclination loop densities. We found that the axis of a 
disclination can be translated to a new position by adding a discrete dislocation line to it. We 
found the basic elastic fields as closed line integral, which confirmed again their nature as state 
quantities. 

The continuous distribution of infinitesimal defect loops was defined by extending Kroupa's 
definition for dislocation loops. The complete correlation to the continuous distribution of defects 
was established by finding the basic plastic fields in terms of the loop densities. It was then simple 
to identify Mura's plastic quantities, as mentioned above, which resolves a difference we had 
with Mura. 

We found that the three independent definitions for defects, i.e., for continuous distributions, 
discrete lines, and infinitesimal loops, do not lead to a single concept of the disclination. Rather, 
the disclination defined in one formulation contains a certain amount of dislocation in the two other 
formulations. So a certain amount of mixing occurs in going from one formulation to another, but 
the amount can be uniquely determined. Our difference with Anthony for example originates from 
this mixing between continuous defect distributions and discrete defect lines. As another example 
we interpreted the discrete defect line in terms of a continuous distribution of defect loops, which 
clarifies our difference with Mura. 

The discrete dipole line was defined by extending Kroupa's definition of the dislocation dipole. 
We find there is a great similarity between disclination dipoles and dislocation lines. 

Then we showed that the general defect theory reduces to the classical dislocation theory when 
the disclinations vanish. 

We clarified the concept of the "dislocation model" of a discrete defect line, introduced by Li 
and Gilman. We found that it is a distribution of dislocation loops identical to that for the corre- 
sponding discrete defect, but without any distribution of disclination loops. In Mura's terms it is 
given by the "plastic distortion" without "plastic rotation." We found that the dislocation model is a 
dislocation wall which terminates on a discrete dislocation lying along the same line as the corre- 
sponding defect line. The dislocation model has the same total displacement, elastic strain, and stress 
as the corresponding discrete defect line. The "elastic distortion" of the dislocation model is 
Mura's "elastic distortion" of the discrete defect line. The great similarity between the dislocation 
model and its corresponding discrete defect line makes it important to distinguish clearly between 
these two concepts. 

We concluded the paper with the "compensated disclination line," which is a constant distribu- 
tion of disclination loops over a surface. In Mura's terms it is given by the "plastic rotation" without 
"plastic distortion." We found that it is a dislocation wall terminating at a discrete disclination 
line, giving no displacement, strain, and stress. The sum of the compensated disclination line and 
the dislocation model gives the discrete defect line. 

In an appendix we have developed a special notation, adapted from Kunin. It is very helpful 
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for the treatment of discrete defects where generalized functions appear, such as the Dirac delta 
function. 

So we have presented a general theory of stationary defects for a linearly elastic, infinitely 
extended, homogeneous body. In future publications we shall specialize these results to isotropy [24] 
and apply them to straight disclinations [25] and disclination loops [26]. 

The major shortcoming of the present treatment might be the use of linear theory. This means 
that in a real solid the resulting fields close to discrete defects will deviate considerably from our 
formulas, but they will become more realistic the further away we are from a defect. This point 
will be more clearly illustrated in the future publications where we obtain specific results for par- 
ticular geometries. However, without the linear assumption we certainly could not have pushed 
the theory as far as we did. This is the price we paid for a fairly complete analytic treatment. 

Within its limitations the present theory is completely self-consistent. Aside from its possible 
intrinsic usefulness, it can be used as a starting point for generalizations, such as dynamics, 
nonlinear effects, couple-stresses, a finite body, or inhomogeneities. 

1 2. Appendix A. The Divergence Theorem and Stokes 7 Theorem 

The rank (also called order by some authors) of a tensor equals the number of subscripts on 
the tensor. In the following let T be a tensor of any rank, where we have suppressed the subscripts. 
The divergence theorem is formulated as follows 



l T ' d H 



TdS,. (Al) 



where the integrations are restricted to the arbitrary volume V and its bounding surface S, which is 
therefore a closed surface. 

Stokes' theorem can be formulated in two ways. The first one is 



I 



€ iik T. s dSt = 6TdL k , (A2) 

where €ij k is the permutation symbol. By using the identity 

£ijk€klm = 8«8j m ~ &i m 8j! (A3) 

we find the second formulation 

j> e klm TdL k = I (T, m dSi — T,tdS m ). (A4) 

Here the integrations are restricted to the arbitrary surface S and its bounding curve L, which is 
therefore a closed curve. We use the right-hand rule in relating the curve L to the surface S. 

13. Appendix B. Delta Functions on Curves and Surfaces 

This appendix is an adaptation of a treatment by Kunin [27]. Let ${x) be an infinitely differentia- 
ble finite function of x, called a test function. We can define the Dirac delta function 8(x) by 



/: 



> b f 0, if x < a, 

8(x-x')<l>(x)dx=\ <£(*'), iia<x < b, (Bl) 

lo, if6<*\ 



where a and b are arbitrary constants. Since 8(x) is a generalized function, the integrals in this 
and subsequent equations are meaningless in the sense of classical analysis. Instead they are to be 
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regarded as a symbolic notation. This notation has been treated in detail by Gel'fand and Shilov 
[28]. The meaning to be attached to the following integrals then follows from the fact that all 
subsequent results are derived from the definition (Bl). 

The three-dimensional Dirac delta function 8(r) is defined by 

S(r-r') = 8( Xl - x[)8{x 2 - x: 2 )d(x 3 - x' 3 )- (B2) 

Let the test function </>(r) be an infinitely differentiate finite function of position r. Then it follows 
from (Bl-2) that 



[ 8(r-r')ct>(r)dV = 



' </>(r'), if r is in V, ( . 

, 0, it r is not in J 7 , 



where dV= dx\dxidx 3 and Fis an arbitrary volume. 

The Dirac delta functions for a curve L, a surface S, or a volume V are defined by 

8i(L)= j 8(r-r')dL' v (B4) 

8i(S) - f 3(r-r')dS f r (B5) 

8(V) - [ d(r-r')dV' (B6) 

We see that 8i(L) and 8/(S), in addition to being delta functions, are also vectors parallel to the 
curve L and normal to the surface S, respectively. From (B3) we see that 



1, if r is in V, 
0, if r is not in V. 



(B7) 



8(V) = 
We now have the following relations 

j &,(L)4>(r)dV=j i 4>(r)dL h (B8) 

j 8 i (S)<l>(r)dV= J <f>(r)dS t , (B9) 

[ b{V)<f>{r)dV= f 4>(r)dV, (BIO) 

where the integrals in the left-hand sides are over all space, and those on the right-hand sides are 
restricted to the curve I., surface S, and volume V, respectively. To prove the first relation, (B8), 
substitute (B4) in the left-hand side, and we have 

[ [ 8(r-r')dL' i <j ) (r)dV= \ [ b{r-r')<j,{v)dVdL\= f 4(r')dL', 

by interchanging integration and using (B3). The proofs of (B9) and (BIO) follow in a similar way. 
By a slight change of variable we can also write (B8-10) as follows 

96 



j 8 ! (L')<f>(r-r')dV'=j i <t,(r-r')dL' { , (Bll) 

J8 l (S')4,(r-r')dV'=f<l>(r-r')dS' i , B12) 

j 8(V')<l>(r-r')dV'= (^(r-r')dV'. (B13) 

If the curve /. crosses the surface S once, we also have the relation 

f faCr-r'^dS^jf^r, (B14) 

where tj is the unit tangent to the curve L, and rij the unit normal to the surface S, at the point of 
intersection, r LS . To prove this relation, note that the integrals in (B14) contribute only at the 
point r ,s . Therefore we can replace the curve L by a straight line tangent to L at r LS , and the surface 
S by a plane tangent to S at r LS . We shall next calculate the integral in (B14) for the special case 
ftj=(001), i.e., when S is the X\Xi plane. Fory = l, 2 the integrals then vanish in agreement with 
(B14). For ; = 3 we have by (B2) and (Bl) 



I J I 8(x i -x[)8(x2-x , j8(x^-x:)dL:dx l dx 2 = I b(x*-x' n )dL[ 

= I 8 (xn -«■ x' 3 ) 7-^-7 dx'z 

J -x I t-s I 

ti 

1 1> I ' 

Here we have used the fact that for a straight line dL'. = (tijt^)dx'^ and that the range of x' 3 is 
(— oo 9 oo) for £3 positive and (°°, — °°) for t :i negative. This result also agrees with (B14) for this special 
case. Hence (B14) holds in a particular coordinate system. By tensor analysis it is therefore also 
true in a general coordinate system. The specific form of (B14) can be derived from the above 
results by the method of Appendix II of reference [29|. From (B4), (B5), and (B14) we have 

* * 1, if i crosses S positively, 

I 8t(L)dSi= f 8i(S)dLi = \ 0, if L does not cross S, (B15) 

I — 1, if L crosses S negatively. 

It is not difficult to generalize these relations to 

/. /. ((f)(r IS ), if L crosses S positively, 

I 8i(L)(/)(r)dSi- I 8 i (S)(f)(r)dL i =\o, if L does not cross S, (B16) 

[— (^(r 7 ^), if L crosses S negatively. 

The derivative of delta functions is defined by switching the operation to the test function 

f 8j(r-r')<f>(r)dV = - J d(r-r')<t>j(r)dV, (B17) 
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as suggested by classical analysis. Now we can write the derivatives of the delta functions (B4-6) 
as follows 

8,j(L) = | 8,j(r-r')dL'., (B18) 

8ij(S) = | 8j(r-r')dS;, (B19) 

8,i(n=fs, j (r-r')#'. (B20) 

From these relations, it follows that 

j 8ij(L)<f>(r)dV=- f <M*0 <^i, (B21) 

| 8,j(S)^(p)dF=-| 0j(r)dSi, (B22) 

| 0>j (F)</>(r)^=- | *j(r)rfF. (B23) 

For example, to prove (B21), substitute (B18) into the left hand side, and we have by (B17) and (B3) 
j j i 8, j (r-r')dL'.<f>(r)dV= jj 8j(r-r')<j>(r)dVdL'. 

= -j i J8(r-r')<t>, J (r)dVdL'. 

— f *j'(r')<tt|. 

The proofs of (B22) and (B23) follow in a similar way. 

The divergence theorem also holds for delta functions and is expressed as follows 

8AV) = -8i(S), (B24) 

where S is the closed surface which is the boundary of V. This theorem is proved showing the 
following relationship: 

J8AV)<t>(r)dV=-j v <t>Ar)dV 
= -j>4>(r)dSi 

= -J8i(S)<t>(r)dV, 

which follows from (B23), (Al), and (B9). Since 4>(r) can be chosen arbitrarily, (B24) follows. We 
also conclude from this relation that 

e ijk 8 k ,j(S) = (B25) 

for a closed surface S. 
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or 



Furthermore Stokes 7 theorem also holds for delta functions as follows 

€ijk8i,j(S) = -8 k (L), (B26) 

6 k i m 8 k (L) = 8 m ,,(S) - 8,, m (S) 9 (B27) 

where L is the closed curve which is the boundary of S. This theorem is proved by showing the 
following relationship: 



je Uf cdij(S)(t>(r)dV=-j* i j k <l>j(r)dS i 
= —<b(f)(r)dL k 
= -jd k (L)<f>(r)dV, 



which follows from (B22), (A2), and (B8). Since 0(r) can be chosen arbitrarily, (B26) follows, and 
(B27) follows directly from (B26). We conclude from (B26) that 

8 k , k (L)=0 (B28) 

for a closed curve L. 

A homogeneous function of degree A is defined by the equation 

f(kx) = k k f(x) (B29) 

for any positive k. We wish to show that 8(r) is a homogeneous function of degree (—3) in r. From 
(B3) we have 

jd(r)(t>(r)dV=<t>(0). (B30) 

By a change of variable and (B3) we have 

\d(kr)<t>(r)dV=k-^8(r)ct>(k-'r)dV 

J J (B31) 

= /c- 3 </>(0). 

Therefore, since <£>(r) is arbitrary 

d(kr) = k-*8(r), (B32) 

which shows the contention. 
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